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Abstract 



The asymptotic giant branch (AGB) phase is the penultimate stage of evolution 
for low- and intermediate-mass stars. Nucleosynthesis products transported from the 
helium-fusing shell to the outer, cooler regions form gas molecules and dust grains 
whose chemistry depends on whether oxygen or carbon is more abundant on the 
surface. Radiation pressure causes the oxygen- or carbon-rich dust to flow outward, 
dragging the gas along. Such outflows inject a significant amount of material into 
the interstellar medium (ISM), seeding new star formation. AGB mass loss is thus 
a crucial component of galactic chemical evolution. The Large Magellanic Cloud 
(LMC) is an excellent site for AGB studies. Over 40,000 AGB candidates have been 
identified using photometric data from the Spitzer Space Telescope Surveying The 
Agents of a Galaxy's Evolution (SAGE) mid-infrared (MIR) survey, including about 
35,000 oxygen-rich, 7000 carbon-rich and 1400 "extreme" sources. For the first time, 
SAGE photometry reveals two distinct populations of 0-rich sources in the LMC: 
a faint population that gradually evolves into C-rich stars and a bright, massive 
population that circumvents this evolution, remaining 0-rich. 

This work aims to quantify the mass-loss return from AGB stars to the LMC, a 
rough estimate for which is derived from the amount of MIR dust emission in excess 
of that from starlight. I show that this excess flux is a good proxy for the mass-loss 
rate, and I calculate the total AGB injection rate to be (5.9-13) x 10^"^ M0yr^^. A 
more accurate determination requires detailed dust radiative transfer (RT) modeling. 
For this purpose, I present a grid of C-rich AGB models generated by the RT code 
2DUST, spanning a range of effective temperatures, gravities, dust shell radii and 

ii 



optical depths as well as a baseline set of dust properties obtained by modeling a 
carbon star, data for which was acquired as part of the spectroscopic follow-up to 
SAGE. AGB stars are the best laboratories for dust studies, and the development of 
a model grid will reinforce future research in this field. 
Advisor: Margaret Meixner. 
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Chapter 1 



Introduction 



1.1 On the Asymptotic Giant Branch 

Low- and intermediate-mass stars (0.8 - 8 M©) spend most of their life in the 
core hydrogen-burning phase (the so-called Main Sequence). This is followed by 
exhaustion of the core hydrogen accompanied by shell H-burning (the "Red Giant" 
phase) and a brief core He-burning stage (the "Horizontal Branch"). The star is now 
left with an inert C/0 core surrounded by a He-burning shell and a dormant hydrogen 
shell, and it has entered the asymptotic gia nt branch (AGB). The H e-burning stage 
( "early- AGB" or E-AGB) lasts for ~ 10^ yr f Vassihadis fc WoodlllQQsl ) and eventually 



reignites the H shell. From this point onward, the star's energy output is due to long 
periods of quiescent shell H-burning punctuated by brief, violent flashes of shell He- 
burning ("thermal pulses") that extinguish the H layer by pushing it to outer, cooler 



regions. The star is now on the t 



5 X 10^ yr. 



lermally-pulsing AGB or TP- AGB. Although brief 



Vassiliadis fc Wood 



19931 ). this is a very important chapter in the star's 
life, and is responsible for a variety of amazing phenomena such as thermal pulsations 
on timescales of hundreds of thousands of years, stellar pulsations on scales of up to 
hundreds of days, and the mixing of nuclear-processed material to the outer regions 
of the star. 

AGB star radii are typically a few hundred R©. The low surface gravity makes 
it easy for stellar pulsations to levitate material to cooler regions leading to the 
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formation of gas molecules and dust grains. Once an appreciable density is reached, 
the dust grains efficiently interact with the incident stellar radiation and are driven 
away from the star, dragging gas molecules along in slow (~ 10 km s^^) winds. 
The wind-forming region typically reaches 10-20 stellar radii, and the circumstellar 
envelope (CSE) formed by the wind can extend to ^10^ R©. These outflows produce 
large mass-loss rates ranging from 10~'^-10~^ M0yr^^ (compare this to the current 
mass-loss rate from the sun, lO"^'* MQyr"^). The dust grains in the CSE absorb 
stellar radiation at UV and optical wavelengths and reradiate it in the infrared (IR), 
making AGB stars some of the brightest sources in the IR sky. Increased formation 
of dust results in progressively more obscuration of starlight until eventually the star 
cannot be detected at visible wavelengths. As the star evolves along the TP- AGB, 
products of shell- burning are deposited onto the inert C/0 core, increasing its mass. 
If the core mass exceeds the Chandrasekhar limit, the star must explode in the form 
of a supernova type 1.5. Since these have never been observed, it follows that stars 
with initial masses greater than the Chandrasekhar limit must lose mass at rates 
faster than the nuclear reaction rates. Calculati ons of the mass-loss rates from the 
brightest AGB stars have shown this to be true (Ivan Loon et al.lll999h . 

The carbon-rich nuclear processed material is mixed into the outer layers of the 
star during the third dredge-up process, which increases the relative abundance of 
carbon compared to the oxygen that is already present in the AGB atmosphere. The 
less abundant element of the two is locked up in the extremely stable CO molecule, 
and the more abundant species decides the chemistry of the CSE. Based on the C/0 
ratio, we can have 0-rich AGB stars (M stars, silicate dust), C-rich AGB stars (C 
stars, carbonaceous dust) and S stars, where the C/0 ratio is close to unity. As 
the star evolves along the AGB, its C/0 ratio gradually increases, goes from an 
M-star to an S-star and flnally ends up as a C-star. In stars more massive than 
about 4 M0(the value is metallicity dependent), the temperature at the bottom of 
the convective envelope is high enough to convert the carbon formed in the CNO 
cycle into N and O, destroying it before it reaches the atmosphere, thus causing it 
to never become a C-rich AG B star. This process is called Hot-Bottom Burning 



( iBoothroyd fc Sackmann 



19921 ). 
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AGB mass loss injects nuclear processed material into the surrounding interstel- 
lar medium (ISM), thus enriching the environment. AGB stars thus directly effect 
the chemical evolution of galaxies as well as contributing to their integrat ed spec- 
tra. Since the mass-loss rate effectively determines the lifetime on the AGB (jWillson 



20001), it is important for understanding stellar evolution along the AGB. Mass loss 



also has important consequences for the luminosity function (LF) e.g. of carbon stars. 
LF information derived from optical data would miss the most luminous (and most 
obscured in the visual) stars. The contribution of these missing stars could be de- 
termined by measuring the mass-loss rates for a large number of stars with a wide 
range of luminosities. The rate of dust production depends on the metallicity of the 
progenitor AGB stars. It is important to measure the mass-loss rates of the gas and 
dust and therefore the gas-to-dust ratio in environments of varied metallicites in order 
to understand AGB mass loss and quantify the rate of mass-loss return to the ISM. 



1.2 The SAGE survey 



The Large Ma gellan i c Clo ud (LMC) is one of our nearest neighbor s, at a distance 



of ab out 50 kpc (iFeastI Il999l ) and a favorable viewing angle of 35° (INikolaev et al. 



20041 ). These factors permit a detailed view of the resolved stellar populations and 
ISM clouds, making it a great laboratory for studying the life cycle of baryonic mat- 
ter. Studies of these populations in the LMC are not hampered by crowding along 
the line of sight or high extinction due to ISM, unlike the Milky Way or the Small 
Magellanic Cloud. The sub-solar metallicity of the LMC is similar to that of the 
ISM during the epoch of peak sta r formation in the Universe (at z~1.5. See, e.g., 



Madau et al. 



1996 



Pei et al. 



19991 ) which makes it a proving ground for star 



Bekki fc Chiba 



orma- 



20051 ). 



tion and galaxy evolution models in the early Universe (see 
In the recent past, the LMC has been the subject of three rn ajor IR surveys. The 
Infrared Astronomical Satellite {IRAS, iNeugebauer et al.lll984l ) s urveyed the LMC in 



four 



Dands (12, 25, 60 and 100 /im), identifying 18 23 point sources (Schwering fc Israel 



1990l ). The Midcourse Space Experiment {MSX, 



Price et al. 



20011 ) imaged the LMC 
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in one band only (8 /im) and identified 1806 poi nt sources (jEgan et al.ll200ll ). More 
recentl y, the 2 micron All Sky Su rvey (2MASS. ISkrutskie et al.ll2006l ) surveyed the 



LMC flNikolaev fc Weinberg! l2000h in the J, H and Kg bands (1.24, 1.66 and 2.16 



fim), and when combined with t he I (0.8 /xm), J, H a nd Kg Deep Near- Infrared Sur- 
vey of the Southern Sky {DENIS, Epchtein et al.ll994 ). it produced a catalog of about 



a million point sour ces. The SAGE {Surveying the Agents of a Galaxy's Evolution, 



Meixner et al. 



20061 ) survey improved on these observations with a 7° x 7° image of the 



LMC using the IRAC (fiher bands centered at 3.6, 4.5, 5.8 and 8 /im) a nd MIPS (24, 
70 an d 160 /xm) cameras on board the Spitzer Space Telescope {Spitzer, 



Werner et al. 



20041 ). The SAGE Epoch 1 Catalog contained about 4 million point sources, includ- 



ing ~ 45 000 dust evolved stars. The 8 fim data was sensitive e nough to be a ble to 

koo& \ 



Blum et al. 



detect the stars with mass-loss rates down to ~ 10^*^ MQyr"^. 
presented color-magnitude diagrams for the evolved stars in the survey, identifying 
about 17500 0-rich and 7000 C-rich sources, and 1200 sources classified as "extreme" 
AGB stars based on their brightness and extremely red colors (these correspond to 
the heavily obscured AGB stars). The color-magnitude diagrams for these sources 
can be used to assess their relative importance to the mass loss budget in the LMC, 
followed with detailed model calculations for their mass-loss rates. A detailed spec- 
troscopic followup to the SAGE survey has been performed (SAGE-Spec, Kemper et 
al. 2009, in preparation) that will help further constrain the mass-loss rates of these 
sources. Some of the bright LMC AGB stars have already been analyzed (see, e.g.. 



van Loon et al. 



1999 



200 



1 



Zijlstra et al.l 120061 ). 



1.3 This thesis 

The main aim of this thesis is to quantify the mass-loss return from AGB stars to 
the LMC. An accurate measurement of the mass-loss return requires detailed radiative 
transfer modeling of the circumstellar dust shell around each AGB star. This can be 
achieved, for instance, by developing a grid of AGB star models that samples the 
relevant range of stellar and dust shell parameters. Towards this end, I present a 
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model grid for carbon stars in this thesis. A companion study by Sargent et al. (in 
preparation) addresses oxygen-rich AGB stars. The thesis is organized as follows: 
in Chapter [2], I calculate infrared excesses for the AGB star candidates identified 
from the SAGE survey, and present empirical relations showing that these mid-IR 
excesses can be used as a proxy for the mass-loss rate. Furthermore, I derive a rough 
estimate for the dust injection rate from AGB stars into the LMC ISM. In Chapter 
[31 I describe a radiative transfer model for the circumstellar dust around OGLE 
051306.52—690946.4, a variable carbon star observed in the SAGE and SAGE-Spec 
studies, using the 2Dust code. This is done in order to derive a baseline set of dust 
properties for use in the carbon star model grid. In this chapter a simple treatment 
for the circumstellar C2ll2 feature is also provided. Chapter H] describes the details of 
a grid of carbon star models spanning a wide range in the relevant parameter space. I 
provide synthetic photometry from these models in the 2MASS and Spitzer bands as 
well as the AKARI and WISE passbands, and present color-color and color-magnitude 
diagrams in order to compare the models with SAGE observations. I also revisit the 
model for OGLE 051306.52—690946.4 in order to test the applicability of the model 
grid. In the future, this model grid will be fit to the entire SAGE AGB candidate 
list in order to derive mass-loss rates for the entire population, therefore enabling 
a better determination of the AGB dust injection rate. Data from the AKARI and 
WISE missions will help tailor our models further towards simulating more realistic 
AGB dust shells, and we will be able to provide a general purpose fitter for galaxy- 
wide data sets of AGB stars. 
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Chapter 2 

The Infrared Excesses of LMC 
AGB Stars 



Note: The results from this chapter are summarized in ISrinivasan et al.l (120091 ) 

2.1 Introduction 



AGB mass loss is believed to be a two-step process: pulsations first levitate material 
above the photosphere, where the cool temperatures result in the formation of dust 
grains. Radiation pressure then drives the dust grains (which are coUisionally cou- 



pled with the gas) outward in an efficient stellar wind (IWickramasing^ 



Goldreich fc Scoville 



1976 



Bowen fc Willson 



1991 



Wachter et al. 



le et al 



1966 



2002I ). An increase 



in luminosity (and hence radiation pressure) must therefore, in general, be accom- 
panied by an increased mass- loss rate (MLR). T his inference i s supported by vari- 



1993; 


Blocker 


1995 


; ] 


3owen & Willson 


1991 


van Loon et al. 


200. 


3)- 



Reimers 



van Loon et al. 



1975 



1999 



Vassiliadis & Wood 


Wachter et al. 


2002; 



The formation of a circumstellar envelope (CSE) is a direct result of the AGB 
mass-loss process. The flux from the CS dust shell appears as mid-infrared (MIR) 
emission in excess of that due to the central star alone. This MIR excess is directly 
related to the rate of mass loss, and is therefore expected to increase with increasing 
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luminosity of the central star. Various studies in the past h ave demonstrated the 
relationship between the MIR excess and the rate of mass loss. 



Skinner fc Whitmore 



(119881 ) showed that the MLR derived from CO rotational transitions was proportional 



Jur 



2 



to the strength of the 9.5 /im silicate dust feature in optically thin 0-rich stars 
(119871 ) derived a relation betwee n the MLR an d 60 /xm excess for O- and C-rich AGB 
stars in the solar neighborhood. iKnapp et al.l (119921 ) found that the 12 nm emission 
from elliptical galaxies was proportional to the measured MLR. lAthey et al.l (120021 ) 
compared the MIR excess emission from 9 galaxies to that of Galactic and LMC AGB 
stars and derived a proportionality relation between the MIR excess and MLR. 

While the study of Galactic AGB stars is hindered by obscuring Galactic dust 
along the line of sight and high uncertainty in distance estimates, the Large Magel- 
lanic Cloud (LMC) is ideal for such observati ons owing to its proxiniity (^ 50 kpc 



Feast 1999) and favorable viewine an 


Kle (35°, 


van der Marel & Cioni 


2001 


) . Early 


LMC surveys ( 


Westerlund et al. 


1978, 


1981 


Rebeirot et al. 


1983; 


Blanco et al. 1980; 


Blanco & McCarthv 


1983; 


Frogel & BlancJ 


199C 


) looked for AGB stars at optical and 



near-infrared (NIR) wavelengths, and hence preferentially detected sources with opti- 



cally thin CS dust shells. The Infrared Astronomy Satellite (IRAS: iNeugebauer et al. 



19841 ) survey of the 



NIR confirmation (Reid et al. 



MC (jSchwering fc Israeli! 1990l ) in conjunction with ground-based 



losing evolved stars (jLoup et al.l 



19901) helped identify seve r al of the brightest mass- 

(119991 ) per- 



1997 



Zijlstra et al.lll99fih . 



Trams et al. 



formed follow-up mid-infrared (MIR) p hotometry and spec troscopy of 57 sources using 



the Infrared Space Observatory (ISO; iKessler et al. 
chemical composition of the CSEs, and 



1996) in o rder to determine the 



van Loon et al. 



(119991 ) used radiative trans- 
fer modeling to derive luminosities and MLRs for these sources. More recently, the 
LMC ha s been surveyec 
(MCPS; Izaritskv et al. 



in th e optical by the Magellanic C 



1997h and in the NIR by DENIS (lEpchtein et al.lll994l ) and 



2 MASS (ISkrutskie et all 
in 



■ouds Photometric Su rvey 



20061). The LMC data f rom t hese two surveys can be found 



Cioni et al.l (l2006l ) and iNikolaev fc Weinberg (l2000l ) , respective- 



Price et al. 



A mid-infrared 



2001). which was 



survey using the Midcourse Space Experiment (MSX; 
four times more sensitive than IRAS, was performed by lEgan et al.l (120011 ). 

As part of the SAGE survey {Surveying the Agents of a Galaxy's Evolution; 
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Meixner et al. 



20061). a ~7° X 7° area of the LMC was imaged in the Spitzer Space 



Telescope fgpifeer I Werner et aLll2004f ) IRAC (3.6, 4.5, 5.8 and 8.0 /xm) and MIPS (24, 
70 and 160 fim) bands. One of the main goals of SAGE was to detect all the evolved 
stars with MLRs > lO"*^ MQyr^^, to characterize the total rate at which material 
is returned to the ISM by dusty evolved stars, and to better understand the physics 
governing mass loss among evolved stars in the LMC. When complete, the SAGE data 
will be ~1000 times more sensitive than MSX, and they will allow a detailed quan- 
titative derivation of the global mass-loss budget from all stellar populations when 
combined with existing and future MIR spectroscopic observations of evolved stars 



[see, e.a 



gun by 



van Loon et al. 



Meixner et al. 



(1200 



19991 . 



2005 



Ziilstra et al.ll2006r) . Continuing the analysis be- 



Blum et al. 



(l2006l ) identified about 32,000 color-selected 
evolved stars brighter than the tip of the red giant branch (11.85 mag in the IRAC 
3.6 fim band), including 17,500 oxygen-rich (0-rich), 7000 carbon-rich (C-rich), and 
1200 "extreme" AGB stars, and presented color-magnitude diagrams (CMDs) of the 
SAGE epoch 1 data. 

Our aim is to investigate the MLRs of the AGB candidates selected from the 
SAGE survey, which will require detailed radiative transfer modeling of the CSE 
around each such star. In order to simplify this effort to model ~ 10*^ stars, we would 
like to constrain the range of some of the input model parameters, such as the total 
source luminosity, the typical temperature of the CS dust and the optical depth at 
a given wavelength. As a first step towards this modeling goal, therefore, in this 
chapter we take an empirical approach and calculate MIR excess fiuxes due to dust 
from these stars. The distance independence of the M IR colors of AGB sourc es makes 
them a good choice for tracing the mass loss (see, e.q.. lWhitelock et al.ll2006l ). but our 
study based on the MIR excesses for sources in the LMC is aided by the fact that the 
distance to the stars in the LMC is essentially the same. Data from the SAGE survey 
also provide an unprecedented opportunity to include a large number of AGB star 
candidates. We will show that the MIR excess can be used as a proxy for the MLR. 
In a very statistical way, we study the overall trends of excess and derived quantities 
such as color temperature of the dust and MIR optical depth with source luminosity. 
In subsequent papers, we will also be presenting the luminosity function and detailed 
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radiative transfer calculations for the AGB star candidates selected in this paper. 
The chapter is organized as follows: in §2.21 we describe the SAGE database and 
our observational sample of AGB stars. Our procedure for calculating IR excesses is 
explained in §2.31 We present our results in §2.4[ and in §2.51 we compare our results 
to previous work and discuss their implications for future AGB studies in the LMC. 



2.2 Data 



The SAGE catalog and archive point source lists from both epochs of observa- 
tions have been delivered to the Spitzer Science Center (SSC), and are available 
for download a t the SSCU . The SAGE epoch 1 point source catalog is discussed by 



Meixner et al. 



(120061 ) and 



Blum et al. 



(120061 ). In this study, we select sources from 
the IRAC epoch 1 archive instead of the catalog. The archive accepts fainter sources 
than the catalog. Faint limits for both epochs are 18.5, 17.5, 15, and 14.5 mag for 
IRAC 3.6, 4.5, 5.8, and 8.0 fim respectively, compared to 18, 17, 15, and 14 mag for 
the catalog. A source is excluded {i.e., culled) from the archive if it has neighbors 
within a 0.5" radius, whereas this radius is 2" for the catalog. This culling procedure 
ens ures the creation of r eliable lists of point sources at the expense of completeness 



see 



Meixner et al. 



20061 ). In order for two sources to appear blended, they both have 



to be detected, and they have to be bright enough. The probability of such blending in 
the 24 fim band is very low, since the 24 fim catalog has very few sources (~ 40, 000) 
compared to the IRAC catalog. However, we do find a few (~40 out of ~ 10^) IRAC 
sources that are matched to more than one MIPS 24 fim source which we corrected in 
our results and in the main catalog. In addition to this culling, a flux value for each 
archive source in any of the IRAC bands is non-null only if its signal-to-noise ratio 
(S/N) is greater than 5. For the catalog, the S/N is 6 for the IRAC [3.6], [4.5], and 
[5.8] fim bands, and 10 for the [8.0] fim band. As a result of these criteria, the archive 
has slightly more sources and more flux values than the catalog, allowing the inclusion 
of more faint AGB candidates. As of the second delivery, version S13, the archive has 



^http://ssc. spitzer. caltech.edu/legacy/sagehistory.html 
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~4.5 million sources compared to ~4.3 million in the catalog. See the SAGE data 
delivery documento for details of the source selection in the catalog and archive. The 
fact that the fluxes of these fainter AGB candidates are more uncertain when using 
the archive list is mitigated by our requirement of 2MASS detections thus affirming 
that the point source is real. The AGB source lists are extracted from a universal 
table which lists the nearest neighbor in the MCPS and IRAC epoch 1 archive to each 
source in the MIPS24 epoch 1 full catalog (~40,000 sources). The matching radius 
is 3", but the nearest neighbor is used in the match. The IRAC archive sources are 
also bandmerged with the 2MASS catalog. 



We classify sources as low- or moderately-obscured 0-rich/C-rich AGB candidates 



based on their location in the Kg versus J-Kg CMEu (jCioni et al 



2006 



Blum et al. 



20061 ). but we exclude stars without H-band fluxes from our list, as our procedure 
for calculating the excesses (see §2.3p relies on an H-band detection. This criterion 
probably excludes the faintest 0-rich AGB stars, and thus does not significantly 
affect our results. As a separation into 0-rich and C-rich che mistries based on near- 



IR col ors is not possible for the heavily obscured f J-[3.6]>3.1. iNikolaev fc Weinberg 



Blum et al 



(2000)) "extreme" AGB candidates, we follow a procedure similar to 
( I2OO6I ) and select these sources based on their location in the [3.6] versus J-[3.6] 
CMD or (when no 2MASS counterpart exists) the [8.0] versus [3.6]-[8.0] CMD. A 
simple trapezoidal integration of the optical U through MIPS24 fluxes is performed 
to estimate the luminosities for all our sources. 
The 



Kastner et al. 



(120081 ) list of objects with spectro scopic classiflc a tions contains 



Kastner et al 



(120081 ) classify 



14 sources classifled in our study as 0-rich AGB stars. 
10 of these as red supergiants (RSGs). We also flnd that we have classifled 21 of 
their sources as extreme AGB stars - most of these have been identifled as possible 
HII regions, none of them are RSGs. Since these are point sources, we interpret HII 
regions as compact HII regions or massive young stellar objects (YSOs). We thus 



^The file SAGEDataDescription_Delivery2.pdf is available at 
http://data.spitzer.caltecli.edu/popular/sage/20080204^nhanced/documents/ 

•^T he UBVI and JHKg magnitudes are dereddened using the same procedure as in ICioni et al 
(I2OO6I) 
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realize that our sample of AGB stars may be contaminated by non-AGB objects such 
as RSGs and YSOs. We estimate this contamination using simple color-magnitude 
cuts. 

The figures in iBlum et al. 



(120061 ) also show that the AGB stars and RSGs are 
not well separated in the MIR CMDs. It is possible that some of the most luminous 
sources in our list are RSGs. We find 556 objects (124 classified as 0-rich, 17 as 
C-rich, and 415 as extreme AGB candidates) more luminous than Mboi = —7.1, the 
classical AGB luminosity limit. It is not unlikely that some of these sources are AGB 
stars. Luminosities above the classical limit can be achieved by AGB stars at the peak 
of their pulsation cycle (as in the case of OH/IR sta rs) or by stars undergoing hot 



bottom burning (HBB, iBoothroyd fc Sackman: 



19921 ). Nevertheless, these numbers 



provide a very conservative upper limit for the RSG contamination in our sample. 

The re ddest sources in our l ist fall in the region of the MIR CMD also populated 
by YSOs. Iwhitnev et all feoOST ) isolate regions in the [8.0] vs [8.0]-[24] CMD occupied 
more densely by YSO models (Figure 3 in their paper). This separation includes a 
stringent cut at [8.0]-[24]~2.2 (corresponding approximately to a 24 fim fiux > the 
8 fim fiux) to exclude AGB stars. We find that 5 0-rich and 24 extreme AGB star 
candidates in our list are part of the IWhitney et al.l (120081 ) list of high-probability 
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Table 2.1. Source List of Color- Selected AGB Star Candidates 



Identifier^ Type'' magU'^ 5magU magB (5magB magV (5magV magi (5magl magj (5magJ magH (5magH magK (5magK 



SSTISAGEIA 


J054938, 


.72- 


-683458. 


,2 


O-rich 


99, 


.99 


99. 


,99 


19, 


.46 


0, 


.04 


17, 


.55 


0, 


.05 


13. 


,76 


0. 


,07 


12, 


.01 


0, 


.03 


11, 


.11 


0, 


.03 


10, 


.78 


0.02 


SSTISAGEIA 


J055530, 


.35- 


-684647. 


,5 


O-rich 


20, 


.46 


0, 


.19 


18, 


.11 


0. 


.03 


16, 


.18 


0, 


.03 


13. 


,95 


0. 


,04 


12, 


.59 


0, 


.02 


11, 


.72 


0, 


.02 


11. 


.48 


0.02 


SSTISAGEIA 


J055420, 


.11- 


-680449. 


,5 


O-rich 


99, 


.99 


99. 


,99 


18, 


.87 


0, 


.04 


16, 


.59 


0, 


.05 


13. 


,53 


0. 


,03 


11, 


.91 


0, 


.02 


11, 


.02 


0, 


.02 


10. 


.71 


0.03 


SSTISAGEIA 


J055729, 


.20- 


-684444. 


,2 


O-rich 


17, 


.97 


0, 


.06 


17, 


.47 


0. 


.04 


16, 


.32 


0, 


.08 


16. 


,38 


0. 


,30 


11, 


.93 


0, 


.02 


11, 


.03 


0, 


.02 


10, 


.75 


0.02 


SSTISAGEIA 


J055321, 


.17- 


-683114. 


,7 


O-rich 


21, 


.04 


0, 


.24 


18, 


.49 


0. 


.09 


16, 


.51 


0, 


.05 


13. 


,87 


0. 


,03 


12, 


.52 


0, 


.02 


11, 


.60 


0, 


.02 


11, 


.37 


0.02 


SSTISAGEIA 


J054522, 


.57- 


-684244. 


,5 


C-rich 


99, 


.99 


99. 


,99 


20, 


.68 


0, 


.06 


16, 


.45 


0, 


.03 


13. 


,82 


0. 


,04 


12, 


.49 


0, 


.02 


11, 


.17 


0, 


.02 


10, 


.38 


0.02 


SSTISAGEIA 


J055650, 


.80- 


-675030. 


,5 


C-rich 


99, 


.99 


99. 


,99 


20, 


.42 


0, 


.07 


16, 


.66 


0, 


.11 


13. 


,56 


0. 


,04 


12, 


.12 


0, 


.02 


11, 


.01 


0, 


.03 


10. 


.35 


0.02 


SSTISAGEIA 


J055835, 


.31- 


-682009. 


,7 


C-rich 


99, 


.99 


99. 


,99 


20, 


.13 


0, 


.05 


16, 


.28 


0, 


.03 


13. 


,55 


0. 


,04 


12, 


.30 


0, 


.02 


11, 


.02 


0, 


.03 


10. 


.20 


0.02 


SSTISAGEIA 


J055311, 


.71- 


-684720. 


,9 


C-rich 


99, 


.99 


99. 


,99 


21, 


.91 


0, 


.17 


17, 


.93 


0, 


.04 


14. 


,56 


0. 


,04 


12, 


.03 


0, 


.02 


10, 


.91 


0, 


.03 


10. 


.14 


0.02 


SSTISAGEIA 


J055036, 


.67- 


-682852. 


.3 


C-rich 


21, 


.61 


0, 


.44 


19, 


.31 


0. 


.05 


16, 


.68 


0, 


.03 


14. 


,42 


0. 


,04 


12, 


.87 


0, 


.03 


11, 


.92 


0, 


.03 


11, 


.52 


0.02 


SSTISAGEIA 


J052742, 


.48- 


-695251. 


,5 


Extreme 


19, 


.08 


0, 


.09 


18, 


.58 


0, 


.09 


18, 


.72 


0, 


.11 


14. 


,86 


0. 


,07 


14, 


.13 


0, 


.07 


12, 


.41 


0, 


.06 


11, 


.08 


0.04 


SSTISAGEIA 


J052714, 


.19- 


-695524. 


,3 


Extreme 


20, 


.04 


0, 


.12 


19, 


.65 


0, 


.07 


18, 


.36 


0, 


.05 


14. 


,45 


0. 


,05 


13, 


.74 


0, 


.03 


12, 


.13 


0, 


.03 


10, 


.85 


0.02 


SSTISAGEIA 


J053239, 


.06- 


-700157. 


,5 


Extreme 


99, 


.99 


99. 


,99 


20, 


.29 


0, 


.20 


17, 


.30 


0, 


.10 


14. 


,81 


0. 


,04 


13, 


.16 


0, 


.03 


11, 


.73 


0, 


.02 


10, 


.73 


0.02 


SSTISAGEIA 


J052950, 


.52- 


-700000. 


1 


Extreme 


18, 


.80 


0, 


.09 


18, 


.28 


0, 


.09 


17, 


.40 


0, 


.14 


16. 


,08 


0. 


,05 


13, 


.96 


0, 


.04 


12, 


.30 


0, 


.04 


10, 


.80 


0.03 


SSTISAGEIA 


J053441, 


.38- 


-692630. 


,7 


Extreme 


99, 


.99 


99. 


,99 


22, 


.78 


0, 


.45 


20, 


.26 


0, 


.11 


14. 


,59 


0. 


,04 


12, 


.36 


0, 


.02 


10, 


.92 


0, 


.02 


9, 


.87 


0.03 



'The electronic version of this table contains all the AGB candidates considered in this study. 
^SAGE Epoch 1 Archive (SAGEIA) designation, including position coordinates of IRAC source. 

''Sources arc classified as O-rich or C-rich based on their J-Kgcolors, or as "Extreme" based on their 2MASS and IRAC colors. For more details, see i|2.2l 
'^Magnitudes and errors in the UBVI, JHKg, IRAC and MIPS 24 fim bands. A value of 99.99 in any band represents either a saturation or a non-detection. 



Tabic 2.1. (Continued) 



mag36 5mag36 mag45 i5mag45 mag58 5mag58 mag80 5mag80 mag24 (5mag24 



10.53 


0.03 


10.72 


0.03 


10.50 


0.04 


10.44 


0.05 


9.86 


0.07 


11.34 


0.03 


11.40 


0.03 


11.28 


0.04 


11.19 


0.04 


10.59 


0.14 


10.47 


0.04 


10.64 


0.02 


10.46 


0.05 


10.40 


0.04 


9.97 


0.07 


10.47 


0.03 


10.60 


0.02 


10.39 


0.03 


10.30 


0.04 


10.01 


0.08 


11.07 


0.07 


11.02 


0.04 


10.90 


0.03 


10.76 


0.03 


10.10 


0.08 


9.67 


0.03 


9.70 


0.03 


9.68 


0.04 


9.30 


0.04 


8.89 


0.04 


9.76 


0.04 


9.86 


0.03 


9.81 


0.04 


9.30 


0.04 


9.26 


0.03 


9.34 


0.05 


9.03 


0.02 


8.80 


0.03 


8.68 


0.03 


8.61 


0.03 


9.51 


0.05 


9.43 


0.04 


9.27 


0.04 


8.99 


0.03 


8.66 


0.03 


11.06 


0.04 


11.08 


0.04 


10.89 


0.04 


10.68 


0.05 


10.21 


0.09 


9.02 


0.05 


8.65 


0.04 


8.39 


0.03 


8.00 


0.03 


7.49 


0.02 


9.62 


0.04 


9.12 


0.03 


8.71 


0.03 


8.24 


0.03 


7.61 


0.03 


9.38 


0.04 


9.33 


0.03 


9.19 


0.04 


8.83 


0.03 


8.46 


0.03 


8.94 


0.04 


8.23 


0.04 


7.55 


0.03 


6.85 


0.02 


6.03 


0.02 


9.06 


0.05 


8.33 


0.05 


7.81 


0.02 


7.29 


0.03 


6.74 


0.02 



YSOs. A more conservative estimate is obtained by looking for sources fainter 
than [8.0]=7 with F2i^rn/ F^^im > 1, which puts the YSO contamination in our hsts at 
about a hundred sources (61 0-rich, 14 C-rich, and 34 extreme AGB candidates). 




1.0 1.5 2.0 2.5 3.0 

J-Ks 



Figure 2.1 A Hess diagram showing the locus of AGB stars in the Kg versus J-Kg 
CMD. Stars brighter than the tip of the RGB (Kg =12, dashed line) are selected as 
AGB c andidates, and the se are classified as 0-rich and C-rich based on the CMD 



cuts in ICioni et al.l (120061 ). The 0-rich sources are bounded by the lines Ki and K2 
(shown as solid line in figure) in their paper, while sources redward of K2 are classified 
as C-rich. The extreme AGB candidates (not shown in this figure) are classified based 
on their 2MASS and IRAC colors. 



Figure [2]T] shows the Kg versus J-Kg Hess diagram for the region popula ted by AGB 
stars. The IRAC and MIPS24 CMDs are shown in Figures O and 



Blum et al. 



( I2OO6I ) noted the presence of a fainter, redder population of 0-rich sources (finger 
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"F" in their Figure 6). The same population can be seen in our IRAC-MIPS24 
CMD (Figure [23]) at [8.0] magnitudes fainter than ~10. We use a magnitude cut at 
10.2, shown in the figure as a sohd fine. Throughout this chapter, we will differentiate 
between the bright and faint 0-rich populations based on this magnitude cut. Almost 
80% of the 0-rich stars in our sample are fainter than [8.0] =10.2. Photometric 
information for a few sources of each type is shown in Table 12.11 The entire list 



of AG B star candidates is available with the electronic version of 

(bood). 



Srinivasan et al. 




Figure 2.2 The three types of AGB candidates (0-rich: blue, C-rich: red. Extreme: 
green) on an [8.0] vs [3.6]-[8.0] CMD. The tip of the RGB is at [3.6]^[8.0]=11.9. 

Figure 12.41 shows the luminosity function for all three classes of AGB candidates. 
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Figure 2.3 The 0-rich (blue), C-rich (red) and extreme AGB (green) stars on the 
[8.0]-[24] versus [24] C MP. The sources b elow the solid line are the fainter, redder 
population mentioned in iBlum et al. fl2006h (The fin ger "F" in their Figure 6). About 



80% of the 0-rich stars in our sample belong to this population. 

The lower luminosity limit {Mf,oi ~ —4 ) for the 0-rich and C-rich AGB candidates 
is set by our color-magnitude cuts which exclude sources fainter than the tip of the 
RGB. The C-rich sources have a tighter luminosity distribution than their 0-rich 
counterparts (there are only a handful of C-rich sources brighter than Mhoi = —7.5, 
whereas the 0-rich distribution falls off around Mboi = —8.5). The range in lumi- 
nosities for the extreme AGB candidates is Mhoi = —4.5 to M^oi ~ —10. The extreme 
AGB candidates thus have the highest luminosities in the sample. However, we have 
insufficient information at this point to say anything concrete about the breakdown 
of these sources into 0-rich and C-rich chemistries. The spectroscopic follow-up to 
SAGE will provide some information about the dust chemistry of these extreme AGB 
candidates. These distributions peak at Mboi ~ —5 (0-rich), Mboi ~ —6 (C-rich) and 
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Figure 2.4 The luminosity function for our sample of AGB candidates (blue: 0-rich, 
red: C-rich, green: Extreme, black: total). The vertical dashed line is the classical 
AGB luminosity limit. 



Mfjoi ~ —6.5 (Extreme). The vertical dashed line shows the classical AG B luminosity 
limit. While this limit can be exceeded by deeply embedded AGB stars ( IWood et al. 
19921 ). 0-rich sources with luminosities higher than Mi,oi = —7.8 are more likely to be 
RSGs, while the brightest extreme AGB stars may be massive YSOs. We will discuss 
the astrophysical implications of these luminosity functions in a future paper. 
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2.3 Procedure 



We estimate the MIR excess emission of the low- and moderately-obscured 0-rich 
and C-rich AGB star candidate populations by comparing their observed SEDs to an 
expected SED for the stellar photosphere. The IR excess in each wavelength band is 
calculated by comparing the total flux from the source observed in that band to the 
flux expected from the central star as prescribed by a "best-fit" model photosphere. 
As the emission from CS dust dominates the MIR flux of extreme AGB stars, we set 
the MIR flux equal to the excess in each band for our extreme AGB star candidates. 
In this work, we are interested in describing overall trends in the AGB parameters 
such as IR excess as opposed to detailed models of each source. To this end we 
choose a single model photosphere for each type of AGB star - one that best fits the 
SED shape of AG B stars with little or no du st. We use the plane-parallel C-rich 

mm and the spherical 0-rich PHOENIX 



MARCS r nodels of iGautschv-Loidl et al. 

models of Hauschildt et al. ( IQQol ) to calculate the photospheric AGB star emission. 
The differences between plane-parallel and spherical models are insignificant as far 
as the resulting broad band AGB star photospheric information is concerned. The 
range of model parameters we search for a "best-fit" are as follows: the thirty-two 
solar-metallicity C-rich AGB models have C/0 ratios between 1.1 and 1.8, effective 
temperatures ranging from 2600 K to 3200 K, and surface gravity logg between about 
—0.76 and 0. The ~200 0-rich AGB models of solar mass and solar metallicity have 
effective temperatures ranging from 2000 K to 4700 K, and logg between —0.5 and 
2.5 in steps of 0.5. Synthetic photometry in the optical, near- and mid-infrared was 
obtained from each of these models by convolving their SEDs with the filter response 
curves of the Johnson-Kron-Cousins UBVI Q , 2MASS JHkJ] , and Spitzer IRAC 



"'The MCPS magnitudes were placed on the Johnson-Kron-Cousins UBVI system. The de- 
tector quantum efficiency curve was obtained from the Las Campanas Observatory website, 
http://www.lco.cl/lco/index.html. Filter profiles for the Johnson U , Harris B, V, and Cousins I 
filters were obtained from the references in Table 9 of lFukugit a et al 

^The 2MASS filter relative spectral responses (RSRs) derived by ICohen et all (|2003h were ob- 
tained from the 2MASS All-Sky Data Release Explanatory Supplement on the worldwide web at 
http: / /www. ipac.caltech.edu/2mass / releases / allsky / doc / sec64a.html. 
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1^ and MIPS 24 yum filters. The oxygen-rich models had spectral information in 
the range 0.1-1000 /im and the convolutions were done directly with the models. In 
order to compensate for the insufficient wavelength coverage (~0.5-25 /im) offered by 
the carbon-rich models, the U and B band fluxes were dropped, and the flux in the 
MIPS24 band was extrapolated to 30 fim assuming a Rayleigh- Jeans falloff. 

As we are most interested in obtaining the correct shape of the SED for a best-fit, 
the model SEDs are scaled to the H-band flux of the median SED of the hundred 
bluest sources in V— Kg color. We use the same model to calculate the excesses for 
both the bright and faint 0-rich populations because the bluest 0-rich candidates 
lack detections in the 24 fim band and thus can not be separated into bright and 
faint populations. The model that comes closest to describing the oxygen-rich me- 
dian SED has surface gravity log (7=0 and an effective temperature T=4000 K. The 
corresponding carbon-rich best-fit model parameters (for a C/0 ratio of 1.3) are 
log (^=-0. 43 and T=3200 K. Figure [231 demonstrates the filter-folding for these best- 
fit models (top), and shows the comparison to the median SED of the bluest sources 
(bottom). There are many more 0-rich candidates that are fainter and bluer than 
most of the C-rich candidates, which might explain the considerable difference in the 
best-fit model temperature and radius of the two types of sources. 

To determine the excess flux X^, in a band centered around frequency u due to 
CS emission, the best-fit model flux must first be scaled to the flux of each source at 
some "pivot" wavelength and then the difference between the observed flux and 
the corresponding scaled model flux in that band must be calculated. Fitting the 
SEDs of all our sources to one best-fit model photosphere is a simple first approach 
towards modeling the CSE around each AGB star candidate. The effects of interstellar 
and CS extinction on our sources is minimal in the 2MASS JHKg bands, which also 
contain the wavelength range corresponding to maximum emission from AGB star 



^The IRAC RSRs are plotted in lFazio et al. I (|2004l) . and were obtained from the Spitzer Science 
Center IRAC p ages at http://ssc .spitzer.caltech.edu/irac/spectralresponse.html 

■^The MIPS (|Rieke et al.ll2004 RSRs were obtained from the Spitzer Science Center MIPS pages 
at http:/ /ssc. spitzer. caltech.edu/mips/spectralresponse. html 
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Figure 2.5 The panels on the top show the normahzed SEDs for the best-fit 0-rich 
(left) and C-rich (right) models in light gray. The response curves for the optical 
UVBI, 2MASS JHKs, IRAC and MIPS24 filters are also shown. Folding the trans- 
mission profiles of each filter into the model SED results in a flux in every filter (red 
diamonds). The bottom panels show the median flux in each band for the hundred 
bluest (in V-Kg color) 0-rich (left) and C-rich (right) sources (black squares) plotted 
over the flux of the corresponding best-fit model (red diamonds) scaled to the median 
flux in the H band. The H band magnitudes of these bluest sources are in the range 
11.99 to 11.36 (0-rich) and 11.95 to 10.32 (C-rich). The 0-rich model is a 1 
model with T=4000 K and log^ = 0, while the C-rich model has T=3200 K, a C/0 
ratio of 1.3, and log (7 = —0.43. 



photospheres^ Almost all of our 0-rich sources peak in the H banc|^ 



while over 



two-thirds of the C-rich sources peak in the Kg band. All of the photospheric models 



^The most obscured sources will suffer from CS extinction even in the NIR bands, but choice of 
pivot wavelength is not an issue for these sources, as they are probably members of our extreme 
AGB list. 

^This is partly due to the fact that the opacity of the H~ ion reaches a minimum in the H-band. 
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available to us exhibit H-band SED maxima. For these reasons, we place the pivot 
wavelength in the H band (centered at 1.65 fim). The excess then depends on the 
source flux and model flux F^""^ according to 

(rpmod \ 
(2.1) 

where the subscript H denotes H-band fluxes. This equation may overestimate the 
MIR excesses from the redder C-rich sources with Fk > Fh- On the other hand, 
using the same best-flt model with a Kg -band pivot would produce underestimates 
for the excesses. For the reddest C-rich source in our sample, the difference in 8 /xm 
excess resulting from choice of pivot is ~ 30%. This is to be regarded as an upper 
bound to the error introduced in the excess determination due to the choice of the H- 
band as pivot. While this will not alter the general trends we discuss in this chapter, 
it will affect our numerical results for higher excesses (mass-loss rates). 

The photometric errors associated with the source fluxes were used to estimate 
the error 5X in the calculated excess. An excess measurement in a wavelength band 
centered around frequency v was deemed "reliable" only if its signal-to-noise ratio 
was greater than 3. In other words, 

5X,, 1 

— - < - 2.2 

X, - 3 ^ ^ 

Figure 12.61 illustrates the effect of selecting sources using this "3-0"" criterion. The 
distribution of sources that are rejected based on this criterion is fairly symmetric 
about zero, with a slight asymmetry on the positive excess side, indicating that our 
cut is conservative. A similar distribution of rejected sources is seen for the C-rich 
candidates, but they are substantially fewer in number. In both cases, sources with 
excesses below ~ 0.1 mJy are rejected. 

Very little can be inferred about the chemical composition of the dust shells around 
the extreme sources without spectroscopic conflrmation, although most of these ob- 
jects are probably C-rich. We flnd that 60 o f our extreme AGB candidates are 



identifled as C-rich in the 



Kastner et al 



(120081 ) study, while only 9 are classifled as 
0-rich. Our procedure for calculating excesses relies on a classiflcation into 0-rich 
or C-rich sources, which is not possible for most of the stars in this list. However, the 
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Figure 2.6 Histogram of the 8 /xm excesses for the faintest 0-rich AGB candidates 
(and thus with the lowest S/N ratio excess emission). The sohd hne is for all sources. 
The dotted line shows sources with S/N>3. The dashed line shows sources with 
excesses considered unreliable and hence were excluded from our study. 
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excess emission from their extremely dusty CSEs dominates over the photospheric 
emission in the MIR, so that we can set the excess equal to the MIR flux to a good 
level of approximation. We find ~ 8200 0-rich, 5800 C-rich and 1400 extreme sources 
with reliable 8 /im excesses, and about 4700 0-rich, 4900 C-rich and 1300 extreme 
sources with reliable 24 fim excesses. Table 12.21 shows the excesses c alculated for the 



f ifteen sample sources shown in Table 12.11 The electronic version of 
( 120091 ) contains all the sources with valid 8 fim excesses. 



Srinivasan et al. 



We estimate the temperature of the CS dust and the optical dep th fr om the 8 fim 



and 2 4 /xm excesses in a manner similar to [Thompson et al.l (120061 ) and iDayal et al. 



(119981 ). The continuum dust emission is modeled as a blackbody at temperature T^, 
with optical depth t^, 



X, ocfi,(l-e-"'')5,(A,Trf) 



(2.3) 



where 



fl,y = 71 



IT 



(2.4) 



is the geometrical dilution factor with R^n the distance from the central star corre- 
sponding to maximum emission due to dust at both 8 and 24 /im. For small optical 
depths, the excess is proportional to the opacity k^. If the emissivity of the dust can 
be modeled by a A^^ power law in the relevant range of wavelengths, we have 



X,(xn,X-'B,{X,Ta) 



(2.5) 



Assuming a power-law emissivity in the MIR ignores any effects due to strong 
absorption or emission features from silicate dust for example, which can be significant 
for the more obscured 0-rich sourc es. The pow er-law index will also depend on 
the dust species in general (see, e.g., 



Hofner 



20071 ). The A ^ power-law dependence 



adequately describes the emissivity at MIR wavelengths for the purpose of this chapter 
(the Spitzer 8 and 24 fim bands will only detect the wings of silicate emission unless 
the sources are OH/IR stars.) The color temperature is calculated by constructing 
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Tabic 2.2. 



Mid-infrared Excess Fluxes 



Identifier 


Type 


X36* 


5X36'' X45 


(5X45 


X58 


(5X58 


X80 


(5X80 X24 


(5X24 






(mJy) 


(mJy) (m.ly) 


(m.lyj 


(m.ly) 


(mJy) 


(m.ly) 


(m.ly) (m.ly) 


(m.ly) 



SSTISAGEIA 


J054938 


72-683458 


2 


O— rich 


4 


546 





622 


1 


918 





315 


2 


338 





287 


1 


335 





210 





457 





051 


SSTISAGEIA 


J055530 


35-684647 


5 


O — rich 





843 





287 





727 





171 





726 





135 





450 





091 





213 





053 


SSTISAGEIA 


J055420 


11-680449 


5 


O-rich 


4 


254 





784 


1 


932 





279 


2 


170 





363 


1 


222 





187 





349 





050 


SSTISAGEIA 


J055729 


20-684444 


2 


O-rich 


4 


351 





531 


2 


413 





292 


2 


740 





230 


1 


693 





190 





326 





050 


SSTISAGEIA 


J055321 


17-683114 


7 


O— rich 


2 


387 





654 


2 


392 





268 


1 


909 





163 


1 


313 





100 





428 





049 


SSTISAGEIA 


J054522 


57-684244 


5 


C-rich 


18 


980 


1 


265 


12 


830 





720 


8 


010 





533 


7 


247 





488 


1 


283 





078 


SSTISAGEIA 


J055650 


80-675030 


5 


C-rich 


13 


170 


1 


375 


7 


974 





645 


5 


065 





509 


6 


534 





504 





606 





048 


SSTISAGEIA 


J055835 


31-682009 


7 


C-rich 


29 


860 


2 


606 


31 


540 





989 


26 


200 





957 


16 


000 





617 


1 


761 





076 


SSTISAGEIA 


J055311 


71-684720 


9 


C-rich 


19 


900 


1 


948 


16 


800 


1 


066 


13 


170 





819 


9 


984 





475 


1 


570 





077 


SSTISAGEIA 


J055036 


67-682852 


3 


C-rich 


1 


080 





494 


1 


209 





257 


1 


308 





207 





953 





160 





238 





051 


SSTISAGEIA 


J052742 


48-695251 


5 


Extreme 


69 


590 


2 


958 


62 


220 


2 


082 


50 


750 


1 


422 


40 


600 


1 


094 


7 


258 





160 


SSTISAGEIA 


J052714 


19-695524 


3 


Extreme 


39 


680 


1 


490 


40 


240 


1 


117 


37 


570 


1 


141 


32 


440 





895 


6 


457 





154 


SSTISAGEIA 


J053239 


06-700157 


5 


Extreme 


49 


920 


1 


958 


33 


430 





812 


24 


320 





834 


18 


790 





459 


2 


950 





068 


SSTISAGEIA 


J052950 


52-700000 


1 


Extreme 


74 


320 


3 


098 


91 


780 


3 


533 


110 


200 


2 


910 


117 


100 


2 


549 


27 


820 





396 


SSTISAGEIA 


J053441 


38-692630 


7 


Extreme 


66 


460 


3 


197 


83 


610 


3 


607 


86 


330 


1 


828 


77 


850 


2 


089 


14 


490 





305 



'^The electronic version of this table contains all the AGB candidates with valid 8 //m excesses. 

''The MIPS 24 /xvci and IRAC 8.0, 5.8, 4.5, and 3.6 fim band excess fluxes and errors in mJy. 
''The errors in the excess fluxes have been calculated by propagating the photometric errors: 



x7 



5F^ 



a look-up table with the excess ratio calculated as per Equation [23] for a wide range of 
temperatures, and selecting the value of temperature from this table that reproduces 
the observed ratio of excesses. Once the color temperature is calculated, Equation [23] 
can be solved for the optical depth. The color temperature and optical depth values 
calculated using these simplified equations are primarily useful in investigating the 
trends of color temperature and optical depth with excess. The equations break down 
for sources with high optical depth [i.e., for extreme AGB stars). The temperature 
derived from a ratio of two broad band fluxes is a simplification that will in turn affect 
the optical depth calculation. We will model the CS shells of our AGB candidates 
considering details such as the wavelength dependence of the emissivity and the choice 
of dust species in a future paper to obtain more precise estimates for the temperatures 
and optical depths. 

Since we are only interested in the overall trend of the optical depth with excess 
fiux, we fix Rin at ten times the stellar radius R^,, which is calculated by scaling the 
radii of the model photospheres to the luminosity of each star. According to models 
of du st condensation, most of the dust has condensed by ~ (see, e.g. 



Hofner 



20071 ). In practice, the dust condensation radii and Rm/R* values will vary, but 



this complication is ignored in this work as we are only interested in a comparative 
study of the MIR color temperatures and optical depths. While the temperature of 
CS dust varies with radius from ~1000 to ~100 K, the (single) color temperature 
calculated here will be dominated by the distance corresponding to the hottest and 
most optically thick dust emission in the CS shell, which is typically close to the 
inner radius of the dust shell. Our values for describe the optical depth of the 
high density inner region of the CS shell. This estimate is insensitive to large but 
cool optical depths. This optical depth will be the dominant component of the line 
of sight optical depth. Table 12.31 shows the color temperatures and optical depths for 
the sources in Tables 12.11 and 12.21 Data for the e ntire set of source s with valid 8 fim 



excesses is available with the electronic version of 



Srinivasan et al 



(120091). 
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2.4 Results 



In this section, we investigate the relation of the MIR excess to other observed pa- 
rameters for the sample including luminosity, color temperature of the CS dust, and 
opacity of the CS dust shell. We focus on the statistical nature of the sample and 
global trends of the results. This approach provides the range of physical parameters 
for the sample and is a first and necessary step to obtain appropriate ranges for model 
parameters for any detailed modeling that may follow for the sample. 

2.4.1 Excess-Luminosity Relations 

The calculated IR excesses show a general increasing trend with luminosity at all 
IRAC wavelengths for O- and C-rich AGB sources. This trend is most obvious at 
8 /im (Figures 12.71 and 12. Sp . While the correlation is similar at the shorter IRAC 
wavelengths, the number of sources with significant reliable excesses is considerably 
smaller. This is expected for the warmer, moderately obscured non-extreme sources. 
The 8 and 24 fim excesses of the extreme AGB candidates also correlate well with 
their luminosities, as is evident from Figure 12. 9[ There is a huge spread in 24 fim 
excesses for the 0-rich and C-rich candidates overall, but an increasing trend is 
apparent for the 0-rich AGB candidates at luminosities above ~ 3 x 10^ Lq. The 
large spread at low luminosities may be caused by significant variations in the MLRs 
during the early stages of AGB stars. The plots of 0-rich sources show the bright 
and faint population sources in light and dark grey respectively. We also display the 
median SED for each population as an inset. For the C-rich and extreme sources, 
we show the median SED for three equally-populated bins in luminosity. Relations 
of the form 

logX,(mJy) = A + B\ogL{LQ) (2.6) 
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Table 2.3. Color temperatures and Optical Depths ^ 



Identifier 




Type 


(K) 


5T 


r(8 fj.m)^ 


<5t(8 ^m) 


t(24 fjm) 


5r(24 fim) 


SSTISAGEIA 


J054938. 


,72-683458.2 


O— rich 


393 


28 


0.813 


0.017 


0.938 


0.003 


SSTISAGEIA 


J055530. 


.00 OO^U^ 1 .0 






1 


u.ouo 


U.UZO 


u.yoo 


u.uuo 


SSTISAGEIA 


J055420, 


11-680449.5 


O-rich 


422 


36 


0.885 


0.012 


0.962 


0.002 


SSTISAGEIA 


J055729. 


,20-684444.2 


O-rich 


502 


48 


0.916 


0.008 


0.972 


0.002 


SSTISAGEIA 


J055321. 


,17-683114.7 


O-rich 


401 


21 


0.743 


0.016 


0.914 


0.003 


SSTISAGEIA 


J054522, 


,57-684244.5 


C-rich 


525 


26 


0.922 


0.004 


0.974 


0.001 


SSTISAGEIA 


J055650. 


,80-675030.5 


C-rich 


814 


84 


0.982 


0.001 


0.994 


0.000 


SSTISAGEIA 


J055835. 


,31-682009.7 


C-rich 


705 


32 


0.953 


0.002 


0.985 


0.000 


SSTISAGEIA 


J055311. 


,71-684720.9 


C-rich 


559 


22 


0.931 


0.002 


0.977 


0.000 


SSTISAGEIA 


J055036, 


,67-682852.3 


C-rich 


446 


53 


0.949 


0.007 


0.983 


0.001 


SSTISAGEIA 


J052742. 


,48-695251.5 


Extreme 


521 


10 


0.830 


0.003 


0.944 


0.001 


SSTISAGEIA 


J052714. 


,19-695524.3 


Extreme 


495 


9 


0.793 


0.004 


0.931 


0.001 


SSTISAGEIA 


J053239. 


,06-700157.5 


Extreme 


559 


11 


0.908 


0.002 


0.970 


0.000 


SSTISAGEIA 


J052950, 


,52-700000.1 


Extreme 


456 


5 


0.601 


0.005 


0.867 


0.001 


SSTISAGEIA 


J053441, 


,38-692630.7 


Extreme 


511 


9 


0.808 


0.003 


0.936 


0.001 



-"^Thc electronic version of tliis table contains all the AGB candidates with valid 8 /im excesses. 

^Color temperatures and related uncertainties, derived from the 8 /xm and 24 /im excesses. The uncertainties are 
calculated by propagating the photometric errors. 

''Optical depths and related uncertainties at 24 fj,m. and 8 /im derived from the color temperature. The uncer- 
tainties are calculated by propagating the photometric errors. 



are fit to the 8 and 24 /xm excesses, resulting in the following best-fit parameters: 



A = -6.7 ± 0.5, B = 1.7± 0.1 (O - rich, 8 /xm) 
A = -9.4 ± 3.2, 5 = 2.1 ± 0.7 (O - rich, 24 /xm, L > 3 x 10^ L©) 



A = -6.6 ± 6.2, B = 1.7± 1.5 (C - rich, 8 /im) 
A = -7.0 ±15, B = 1.6 ± 3.5 (C - rich, 24 /im) 
A = -5.0 ± 0.3, 5 = 1.5 ± 0.1 (Extreme, 8 /xm) 
A = -8.3 ± 0.7, 5 = 2.1 ± 0.2 (Extreme, 24 /xm) 

The poor fit to the 24 /im excess for the C-rich candidates reflects the huge spread 
of 24 /xm excesses for these sources. 



2.4.2 Color Temperature and Opacity 

Figures 12.101 and 12.111 show the variation of the color temperature and 24 /xm optical 
depth with 24 /im excess for all three types of AGB candidates. As expected from 
their redder [8] — [24] color, the members of the faint 0-rich population tend to have 
cooler temperatures and higher optical depths than the bright population, and they 
also show a more pronounced variation with excess. There is a huge spread in the 
C-rich AGB color temperatures. The maximum color temperature is higher than in 
the 0-rich case, suggesting that the carbonaceous dust grains become hotter than 



the silicates. AGB winds composed of oxygen-rich comp ounds are less efficient at a 



b- 



sorbing visible photons than those that are carbon-rich (jWallerstein fc Knapplll998l ). 
Carbonaceous grains are efficient absorbers of optical photons and are highly emissive 
at IR wavelengths in comparison to silicates, thus they are more likely to reach higher 
temperatures. We fit the color temperatures for the extreme AGB candidates with a 
power-law relation of the form 

logT,(K) = A + 51ogX,(mJy) 
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10^ 10* 10^ 10* 

Lumlnoalty (Lg) 



Figure 2.7 The 8 fim (top) and 24 fim (bottom) excesses plotted against luminosity 
for the 0-rich sources. Sources from the bright and faint populations (light gray and 
dark gray circles respectively) are shown. The solid lines show, respectively, power- 
law fits to all the sources with rehable 8 /xm excesses (top) and sources with rehable 
24 iim excesses brighter than 3 x 10^ L© (bottom). Representative error bars are also 
shown. Inset: the median SEDs of the bright and faint populations. 
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Lumlnoalty (Lg) 



Figure 2.8 The 8 fim (top) and 24 fim (bottom) excesses plotted against luminosity 
for the C-rich sources. The solid lines are power-law fits to the excess-luminosity 
relations. (The fit is very poor in the 24 /im case) Representative error bars are 
also shown. Inset: the median SEDs for three equally-populated luminosity bins 
(solid line: least luminous one-third of the sample, thin dashed line: intermediate 
luminosity, thick dashed line: most luminous). 
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LiJ 

1 



10^ 10^ 
Luminosity (Lq) 



10® 




Luminosity (L©) 

Figure 2.9 The 8 fim (top) and 24 fim (bottom) excesses plotted against luminosity 
for the extreme sources. The solid lines are power-law fits to the excess-luminosity 
relations. Representative error bars are also shown. Inset: the median SEDs for 
three equally-populated luminosity bins (solid line: least luminous third of sample, 
thin dashed line: intermediate luminosity, thick dashed line: most luminous). 



31 



T 1 1 I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 I I I I I 




I I I I I I I I I I I I I 

0.1 1 10 100 1000 

24 fun Excess (mjy) 

Figure 2.10 The color temperature derived from the ratio of the 8 fim and 24 /xm 
excesses plotted against the 24 fim excess for all three types of AGB candidates 
(blue: faint 0-rich, cyan: bright 0-rich, red: C-rich, green: Extreme). The cross is 
a representative error bar. 
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Figure 2.11 The 24 /xm optical depth derived from the 8 /im and 24 jjm excesses 
for all three types of AGB candidates. The color coding is the same as for Figure 
12.101 Power-law fits to the optical depth-excess relation are shown (dotted: faint 
0-rich, dashed: bright 0-rich, dot-dashed: C-rich, solid: Extreme). The cross is a 
representative error bar. 
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and obtain the following best-fit values: 



A = 2.9± 0.4, B = -0.19 ± 0.02 (Extreme) 



while a similar relation, 



logr24, 



A + BlogX,(mJy) 



(2.7) 



when fit to the 24 jjm optical depths, gives 



A 



A 



A 



0.49 ± 0.12, B = 1.6± 0.1 (O - rich, faint) 
-1.2 ± 0.3, B = 0.9± 0.2 (O - rich, bright) 
A = -1.7 ± 0.3, B = 1.2± 0.2 (C - rich) 
1 = -1.8 ± 0.1, B = 0.79 ± 0.09 (Extreme) 



Figures 12.121 and 12.131 show the variation of the 24 /im optical depth with luminosity 
and color temperature respectively. There is no trend in the optical depth of the O- 
rich candidates with luminosity. For the C-rich sources, there is no correlation, but 
there appears to be an absence of sources with low optical depth at higher luminosities. 
For the extreme AGB stars, there appears to be a positive correlation of higher 
optical depth with higher luminosities. The optical depth decreases monotonically 
with the color temperature for all three types of AGB candidates. CS dust shells 
have temperature gradients from ~ 1000 K in the interior to ^ 100 K in the outer 
regions. With increasing optical depth, therefore, we are probing the progressively 
cooler, outer regions of the AGB star candidates. 



The 8 and 24 /im band excess fluxes show an increase with increasing luminosity of 
the central star. For each type of AGB candidate, the slope of the excess-luminosity 
relation is similar across the four IRAC and MIPS24 bands, though in the IRAC 8 
/im and MIPS24 bands, excesses are "reliable" only beyond ~0.1 mJy (see Figure 
12. 6p . The increase of excess with luminosity is consistent with the MLR-luminosity 



2.5 Discussion 
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Figure 2.12 The 24 fim optical depth as a function of luminosity for all three types 
of AGB candidates. The color coding is the same as for Figure 12.101 
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Figure 2.13 The 24 /im optical depth plotted against color temperature for all three 
AGB types. The color coding is the same as for Figure 12.101 
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relation found by_ 
to tL (see, e.q. 



van Loon et al. 



Ivezic fc Elitzur 



1999|). The mass loss-rate is roughly proportionate 



19951 ). so we expect that the MLR increases with 



increasing luminosity as long as the optical depth does not decrease faster than ~L^^. 
This suggests that the excess is a good reflection of the MLR. More quantitative 
comparisons of MLR for this whole sample with luminosity is beyond the scope of 
this study, but this will be addressed in a future paper. 

The uncertainties in photometry alone cannot account for the considerable spread 
in the calculated excesses. Stars on the AGB suffer from variability and episodic 
mass loss, and our observations of these stars may be capturing their fluxes at dif- 
ferent epochs in their variability cycle, or during episodes of increased or diminished 
mass loss. At a given luminos ity, these ef f ects i ntroduce a variation in the MLR and 



hence also affect the excess. 



Vijh et al.l (120091 ) show a median variability index of 



about 3 for the 0-rich and C-rich AGB sources and about 5 for the extreme sources. 
This would correspond to a change in the excess fluxes of about the same factor. The 
color temperature depends on the ratio of the excesses, and the effect of variability is 
probably weaker. The effects of grain chemistry, metallicity of the environment, and 
the C/0 ratio of these sources have not been accounted for in this study, as this re- 
quires spectroscopic determination. At a given luminosity, there is also a dependence 
on progenitor mass - a degeneracy between a less evolved, more massive star and a 
more evolved, less massive star arises. Such degeneracy effects increase the spread at 



the lower luminosity end (only the most r aassive stars evolve to high luminosi 



Observations for Galactic AGB stars (IGuandalini et al. 



200 



: ies) . 



Busso et al. 



2003) 



show that C-rich AGB stars are in general more obscured than their 0-rich counter- 
parts, but this is partly due to the higher opacity of carbonaceous dust. Note that 
extreme 0-rich AGB stars {i.e., OH/IR stars) are even more obscured than these 
C-rich sources. In the case of the extreme sources, there is a high optical depth re- 
gardless of chemistry of the CSE. As an AGB star evolves, more mass in the form of 
dust is deposited into the surrounding shell, increasing its opacity and causing further 
reddening of stellar light and stronger emission at longer wavelengths, corresponding 
to cooler color temperatures. The color temperature and optical depth variations for 
extreme AGB stars in Figures [2.101 and 12.111 support these claims. The existence of a 
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redder, fainter population of 0-rich sources first seen in the IRAC CMD is also appar- 
ent in plots of the excess, color temperature and optical depth. The faint population 
is in general cooler and more obscured than the bright population. The SAGE survey 
is able to differentiate between these two types of sources for the first time. The 
bright population corresponds to the young, most massive AGB stars that prevent 
the dredge-up of carbon by undergoing hot-bottom burning (part of sequence G in 
Nikolaev &: Weinberd (120001 ) ) . This population also emerges in the model isochrones 



of iMarigo et al.l (l2008l ) - they find that the thermally-pulsing (TP) phase in their 
isochrone of age log(t/yr) = 8.2 is well-developed and populated by 0-rich stars 
undergoing HBB. The faint population consists of AGB stars in which HBB does 
not occur, resulting in a s mooth transit i on fro m 0-rich to C -rich chemistry. The 
l og(t/ yr) = 9.1 isochrone of iMarigo et al.l (120081 ) (sequence J in lNikolaev fc Weinberg 
(120001 )) represents the typical evolutionary phase corresponding to this population. 

Our cr-clipping method of recovering sources with reliable excesses can be used to 
compare our r esults with the expec ted lower-limit to measurement of MLRs by the 
SAGE survey jMeixner et al.ll2006f ). The mid-infrared CMDs in Figures [2l4l EUSl 
and 12.161 show the sources with reliable 8 /im excesses as black circles. We detect 
0-rich and C-rich sources with reliable excesses up to within ~ 0.1 m ag of the tip of 



Meixner et al. 



the R GB ([8.0]~11.9). which is below the detection limit mentioned in 
(120061) ([8.0]=11.0) for measuring significant mass loss. The extreme AGB stars in 
our sample are all well above this detection limit, showing that the SAGE survey 
detects all the extreme mass-losing AGB sources in the LMC. 

We estimate the current LMC AGB star mass- lo ss budget using a meth od analo- 



gous to 



Blum et al. 



van Loon et al 



(liggj For spec- 



(120071 ). The MLRs obtained by 
troscopically identified AGB stars are plotted against their SAGE 8 /xm excess fluxes 
in Figure 12.171 The figure also shows power-law fits of the form 



logM(10-^ Moyr-^) = A + 51ogX^(mJy) 
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Figure 2.14 [8.0] vs [3.6]-[8.0] and [8.0] vs [8.0]-[24] CMDs for the 0-rich AGB stars 
in our list with 8 and 24 fim detections (gray circles). The sources with "reliable" 
excesses are superimposed (black circles). The tip of the RGB is at [8.0]=11.9 (dashed 
line) 
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Figure 2.15 Same as Figure [2.141 for our C-rich AGB candidates. 
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Figure 2.16 Same as Figure [2 .141 for our extreme AGB candidates. 



to all three types of sources with the following best-fit parameters: 



A 



A 



A 



1.5, B = 1.0 (0-rich) 
•1.7, B = 1.1 (C-rich) 
2.8, B = 1.7 (Extreme) 



These relations are then used to derive MLRs for all the sources in our lists with 
reliable excesses. Figure 12.181 shows the total dust MLR as a function of LMC AGB 
luminosity. The extreme AGB stars, despite their comparatively low numbers, are 
the most significant contributors to the AGB mass-loss budget in the LMC. We find 
the dust injection rates from 0-rich, C-rich and extreme AGB stars to be 0.14, 
0.24, and 2.36 xlO^^ MQyr"^ respectively, which puts the estimate for the total dust 
injection rate from AGB stars at 2.74 x 10"^ MQyr"^. While stars with luminosities 
between Mf,oi = —7.1 and M^oi ~ —7.8 may be highly embedded AGB stars, sources 
brighter than Mhoi > —7.8 are definitely supergiants (Sloan et al., in preparation) 
The luminosity distribution (Figure [2]3]) shows that there are very few sources in this 
range, but it is clear from Figure 12.181 that the contribution of the red supergiants 
to the dust injection rate is comparable to that of the AGB stars. This introduces 
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Figure 2.17 The mass- loss rates from Ivan Loon et al.l (119991 ) for 0-rich (triangles), 
C-rich (filled squares) and extreme (circles) AGB sources plotted against their 8 fim 
excesses. The lines show power-law fits for the dust mass-loss rate-excess relation for 
all three types of AGB candidates (solid line: 0-rich, dashed line: C-rich, dot-dashed 
line: Extreme). 



a significant uncertainty in the value for the AGB dust injection rate derived in this 
section. 

Calculation of the total (gas+dust) injection rate requires knowledge of the gas:dust 
ratio, ip, in the circumstellar shells of AGB stars. The lower metallicity of the LMC 
will result in a higher s ilicate gas: du st ratio but the carbon gas:dust ratio may be sim- 
ilar to Galactic values (jHabing||l996l ). Keeping this in mind, we choose gas:dust ratios 
of 200 and 500 for C -rich and 0-rich sta rs respectively^"!. The C-rich value is similar 
to that obtained by 



Skinner et al. 



(119991 ) for the Galactic carbon star, IRC+10 216. 
The gas MLR from the 0-rich and C-rich AGB candidates is then 7xl0~^ Moyr"^ 
and 4.8x10"'^ MQyr"^ respectively. The chemical identification of the extreme AGB 
sources is not possible from our data, but assuming ip=200 and '?/'=500 for these stars 
will provide lower and upper limits. The gas MLR for the extreme AGB population is 



in the range 4.7-11.8x10"'^ Moyr^^ Thus, the total AGB gas MLR is (5.9-13) x 10 



-3 



-"^"The work of van Loon et al.l ( 20081) suggests that ■0 scales in proportion to the metaUicity for 
both 0-rich and C-rich AGB stars, but this effect can not be accounted for in our simple estimate. 
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Figure 2.18 The cumulative dust mass-loss rate as a function of luminosity for the 
0-rich (thin solid line), C-rich (dashed line) and Extreme (dot-dashed line) AGB 
stars in our lists. The thick solid line is the total AGB dust mass-loss rate as a 
function of luminosity. The extreme AGB stars are the major contributors to the 
mass-loss rate, at 2.36 xlO""^ Mgyr"^. The total dust mass-loss rate is 2.74xl0~^ 
M0yr~^. The vertical dotted line is the classical AGB luminosity limit. As is evident 
from the figure, there are very few, if any, carbon-rich AGB candidates above this 
limit. However, deeply embedded extreme AGB stars and 0-rich stars undergoing 
hot-bottom burning can exceed this luminosity limit. 
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MqJt-^. We will improve on our estimate with our follow-up radi ate transfer mod- 
eling of the dust shells around these stars using the 2DUST code (lUeta fc Meixner 



20031 ). The results of the present study will help constrain the range of optical depths 



and grain temperatures. The optical depth is an important input parameter for the 
code. 

Taking into account the effects of variability on the excesses as discussed previ- 
ously, we estimate that the MLR of an individual source can vary on average by a 
factor of ~3 for the 0-rich and C-rich sources and ~15 for the extreme AGB can- 
didates. Due to the large number of 0-rich and C-rich sources, the errors of the 
flux variations on the cumulative MLR estimates would probably cancel out. For the 
extreme AGB candidates, the larger effects of the variability and the smaller numbers 
may result in a significant change in the cumulative MLR. The magnitude of such 
an effect could be calibrated with an IR variability study of the LMC AGB stars to 
determine their periods and hence phase-correct the fluxes and excesses, but no such 
study exists for our complete sample. 

Our C-rich AGB star gas MLR of 4.8x10"'^ MqJt^^ is comparable to the value, 



6x10 ^ M0yr \ found by 



Matsuura et al. 



(120091 ). However, we note that 



Matsuura et al. 



(120091 ) base their measurement on MLR versus color relations and include the extreme 
C-rich AGB stars, identified by IRS spectroscopy, as well as the color classified C-rich 
AGB stars. Our excess versus MLR relation method is two orders of magnitude more 
sensitive to low e r ML R sources (dust MLRs of ~ 10"^^ MQyr^^) compared to the 



Matsuura et al. 



(120091) method (dust MLRs of ~ lO^^ Moyr^^). It is therefore fortu- 
itous that our estimates are so similar and our exclusion of the highest MLR extreme 
C-rich AGB stars is balanced by the inclusion of more low MLR C-rich AGB stars. 
Furthermore, the total mass injection rate of C-rich mass loss is probably higher than 
both these estimates, but definitive numbers will require a clear identification of what 
fraction of the extreme AGB stars are C-rich. 

Our estimate for the total LMC AGB gas MLR is higher than those obtained 
for two local group d warf irregulars , WLM (jjackson et al.l l2007al . (0.7-2.4) xlO~^ 
Meyr"^) and IC 1613 (I.Tackson et al.ll2007bl . (0.2-1.0) x 10"^ Meyr"^). As their val- 
ues for the MLRs of individual stars are in good agreement with the estimates of 
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individual stars by 



van Loon et al. 



(119991 ) for the LMC, the disparity between our 
results and the values for WLM and IC 1613 could arise due to various reasons. The 
total MLR return due to AGB stars depends on the total number of AGB stars in a 
galaxy. The total stellar mass of the LMC (~ 3 x 10^ M^ 



is two orders o f magni tude greater than those estimated by 



van der Mare 



Jackson et al 



et al 



20021 ) 



Jackson et al. 



(l2007aD and 

J2007bl ) for WLM and IC 1613 (1.1 x 10^ Mr;, and 1.7 x 10^ Mr;, respec - 
tively). Moreover, the age of the LMC bar is 4-6 Gyr ( ISmecker-Hane et al.ll2002l ). 
which is optimal for an enhanced AGB populati on at present times, c ompared to the 



slightly older stellar population of WLM (9 Gyr, I Jackson et al.ll2007al ). Based on the 
differences in stellar masses and ages, we expect that the LMC has considerably more 
AGB stars than either WLM or IC 1316, leading to a higher total MLR from AGB 
stars. However, the differences in the MLR determination may also increase our MLR 
comp ared to theirs. Our IR exces s method is ab l e to det ect sources with small MLR, 



while 



van Loon et al. 



fll999f ) and 



Jackson et al. 



(I2007al ) derive their estimates from 



the highest mass-losing sources alone. In this sense, ours is a more complete es t imate 

Jl999[ ) 



van Loon et al. 



of the total MLR. On the other hand, our interpolation of the 
MLRs may not be appropriate for sources with very low excess. Thus, while the 
excess method is able to identify sources with low MLRs, the derived MLRs may be 
overestimates. 

The curr ent star formation r ate (SFR) of the LMC is estimated to be about 



0.1 MqJt-^ flWhitnev et al 



20081 ). which is an order of magnitude higher than our 



calculated AGB MLR. Thus, the current star formation rate is not sustainable, as- 
suming the LMC is a closed box system, and that AGB stars are the only means 
of replenishing the ISM. However, the most massive stars {e.g., red supergiants, su- 
pernova explosions) may also contribute to the mass-loss return, and these need to 
be accounted for to complete the census of the mass budget return to the LMC. In 
addition, interactions with the Small Magellanic Cloud may cause infall or loss of 
mass that would need to be considered in the mass budget of the ISM. 
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2.6 Summary and conclusions 

We classify evolved stars in the SAGE Epoch 1 Archive based on their 2MASS and 
IRAC colors, and use photospheric models to estimate the infrared excess emission due 
to circumstellar dust. We obtain about 8200 0-rich, 5800 C-rich, and 1400 extreme 
AGB sources with rehable 8 fim excesses (SNR > 3). The corresponding numbers in 
the 24 fim band are about 4700, 4900, and 1300 respectively. The excesses increase 
with an increase in luminosity in all four IRAC bands as well as the M1PS24 band. 
We use the 8 and 24 /xm excess fluxes to derive dust color temperatures and optical 
depths, and we observe that higher excess fluxes correspond to cooler temperatures 
and optically thicker dust shells. These quantities for our list of AGB candidates with 
valid 8 fjLUi excesses are available in the form of electronic tables. We also estimate 
the present day AGB mass-loss budget in the LMC by comparing modeled mass-loss 
rates with our excesses estimates to find that the extreme AGB stars are the most 
significant contributors to mass loss in the LMC. Our data also suggests that the rate 
of dust injection from red supergiants to the LMC ISM is comparable to the total dust 
injection rate from the AGB stars, which we calculate to be about (5.9 — 13) x 10~^ 
Moyr-i. 
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Chapter 3 

Modeling The Circumstellar 
Envelope of OGLE 
051306.52-690946.4 



3.1 Introduction 



In order to develop a grid of LMC carbon star models, it is essential to first identify 
the composition of the dust driving the mass loss in C-rich AGB stars and validate the 
properties of carbonaceous dust produced. In this chapter, I present a 2Dust radia- 
tive transfer model for the circumstellar dust around the variable carbon star OGLE 
051306.52— 69 0946.4. This sour c e has been studied as part of rnany recent LMC 



surveys [e.g., 



2007 



Zebruh et al. 



Meixner et al 



date in 



GroenewegenI (120041 ). while 



2001 



Epchtein et al 



1999 



Cutri et al 



2003 



Kato et al. 



2006). This ob 



ect has been listed as an obscured AGB candi- 



Blum et al. 



fl2006h and 



Srinivasan et al. 



fl2009h 



classified this star as an "extreme" AGB candidate (J-[3.6]>3.1). The SAGE-Spec 
spectrum for this source showed signatures of carbonaceous dust, confirming that 
OGLE 051306.52-690946.4 is an extreme AGB with C-rich chemistry. 

The choice of OGLE 051306.52—690946.4 for this study was motivated by the 
availability of both SAGE photometry and SAGE-Spec spectroscopic data, and the 
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good agreement (;^10%) between the two. The carbon star candidates chosen for the 
SAGE-Spec study were required to be infrared-bright. Thus, OGLE 051306.52—690946.4 
is redder than most carbon stars in the SAGE sample. Despite this fact, the absence 
of pecuhar features in its spectrum allow us to consider it as a good testbed for the 
dust properties around LMC C-rich AGB stars. 

I begin by describing observational data for OGLE 051306.52—690946.4 from vari- 
ous studies in §3.2[ In §3.31 1 detail the method employed to create a radiative transfer 
model for the circumstellar dust as well as a model for the molecular features observed 
in the mid-infrared (MIR) spectrum. I present the results of these models for OGLE 
051306.52—690946.4 in §3.41 and discuss the implications of these results for the dust 
properties to be used in the model grid in §3.51 



3.2 Observations 

Simultaneous measurements in the optical, NIR and MIR wavelengths are essential 
for studying the spectral energy distribution (SED) of a variable source. Although 
SAGE data has been combin ed with photometry from the optical M a.gellanic Clouds 
Phot ometric Survey (MCPS; 



Zaritsky et al. 



19971 ) as well as 2MASS fISkrutskie et al. 



20061 ). these observations were not concurrent with SAGE. In order to constrain the 
optical and near-infrared (NIR) variability, I discuss OGLE and IRSF observations 
for OGLE 051306. 52-690946.4 in In 03X2) and E^l I describe Spttzer Space 

Telescope ( Spitzer AGehrz et al.ll2007l ) IRAC/MIPS photometry and IRS spectroscopy 



for OGLE 051306.52-690946.4 (SAGE J051306.4 0-690946.3) ob t ained as part of the 
SAGE and SAGE-Spec studies, respectively (See iMeixner et al. 
et al. (2009, in preparation) for program overviews). 



(l20o3) and Kemper 



3.2.1 Variability in optical and NIR data 

OGLE 051306.52— 690946.4 was observed as p art of the Optical Gravitational 



Zebruh et al. 



Lensing Experiment (OGLE-II, 
The catalog contains BVI light curves for 68 000 variable objects. 



20011 ) survey of th e Magellanic Clo uds. 



Ita et al. 



(|2004h and 
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GroenewegenI (12004) crossmatched th ese data with the IRSF LMC S urvey ( 



20071), DENIS (lEpchtein et aL 



Kato et al. 



19991) and 2MASS All-Sky Release fICutri et al. 



200^ 



NIR catalogs and fit the li ght curve s to ob tain variability periods and amplitudes 

(361 d) 



The period determined by 



Ita et al. 



(2004) for OGLE 051306.52-690946. 



GroenewegenI (|200J) (358.6±0.136 d). 



Ita et al. 



was in good agreement with that of 
( I2OO4J ) also reported a mean I magnitude of 17.483 mag and an amplitude of 1.865 
mag. 

SAGE data is also matched with 2MASS as well as the IRSF survey of the LMC 
( IKato et al.ll2007l ). which contain JHKs photometry for OGLE 051306.52-690946.4. 
While the NIR information from these surveys is useful for qualitative confirmation, 
the fact that they are not coeval with SAGE observations leads us to expect a slight 
disagreement between the spectral energy distributions (SEDs) in the NIR and MIR. 
We model the source using SAGE Epoch 1 photometry due to its agreement with the 
observed spectrum, and use the optical and NIR data as to constrain the shape of 
the resulting SED. 



3.2.2 SAGE photometry 



The SAGE (Surveying the Agents of a Galaxy's Evolution, iMeixner et al.ll2006l ) 
survey imaged a ~ 7° x 7° area of th e Large Magellan ic Cloud (LMC) with the 
Spitzer Infrared Array Camera (IRAC; iFazio et al.ll2004l^ 3.6. 4.5, 5.8 and 8.0 fim) 

24, 70 and 160 /xm). 



Rieke et al 



2004 



and Multiband Imaging Photometer (MIPS; 
Two epochs of observations separated by three months were obtained to constrain 
source variability. Spectroscopic follow-up observations using the Infrared Spectro- 
graph (IRS) aboard Spitzer were performed as part of the SAGE-Spectroscopy survey 
(Kemper et al. 2009, in preparation) on the brightest objects observed in the SAGE 
survey. 

IRAC and MIPS epoch 1 archive photometry for OGLE 051306.52-690946.4 
is available as pa rt of SAGE data delivered to the Spitzer Science Centei0 (SSC). 

(I2OO6I ) presented an overview of the evolved star population in the SAGE 



Blum et al 



^Data version S13 available on the SSC website, http: / /ssc. spitzer. caltech.edu/legacy/sagehistory.html 
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[3.6]-[B.O] 
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Figure 3.1 The [8.0] vs [3.6]-[8.0] CMD with OGLE 051306.52-690946.4 (filled dia- 
mond) superimposed o n the 0-rich (blu e ), C- rich (red) and extreme (green) AGB 
candidates identified by ISrinivasan et al.l (120091 ) . 



data of the Large Magellanic Cloud (LMC) and identified asymptotic giant branch 
(AGB) star candid ates based on their loc ations in NIR and MIR color-magnitude 
diagrams (CMDs). ISrinivasan et al.l (120091 ) extracted AGB star candidates from the 
SAGE IRAC Epoch 1 Archive list using similar color cuts (see paper for details). Fig. 
O shows the location of OGLE 05130 6.52-69094 6 .4 on a [8.0] vs. [3.6]-[8.0] CMD. 



By comparing data from both epochs, IVijh et al 



(120091 ) identified variable sources 



in the LMC and found that AGB candidates make up ~ 81% of this population. I 
use both epochs of SAGE photometry for OGLE 051306.52—690946.4 for the present 
study, as shown in Fig. 13. 2[ 



3.2.3 SAGE-Spec data 

Basic Calibrated Data products for OGLE 051306.52—690946.4 were obtained 
from the SSC. The spectra were reduced using techniques described by Kemper et 
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al. (2009, in preparation). OGLE 051306.52—690946.4 was observed in both spectral 
orders of the Short-Low (SL; 5.2-14.3 /im, R~60-120) and Long-Low (LL; 14.3-13.7 
/im, R~60-120) instruments on board Spitzer at two nod positions. Data from the 
S15.3 and S17.2 pipehnes for SL and LL, respectively, were obtained from the SSC for 
OGLE 051306.52-690946.4 (AOR # 22415360). SL sky subtraction was performed 
by subtracting the observations in one order from the other while keeping the nod 
position fixed, while LL sky subtraction involved subtracting one nod position from 
the other keeping the order fixed. Bad pixels on the detector arrays were replaced 
with values determined by comparing to the local point-spread functioij^. 

Point source signal extraction used apertures at each wavelength whose width 
increased in proportion to the wavelength; this involved using the profile, ridge and 
extract moduleqj. Flux densities were obtained by calibrating the raw signal extrac- 
tions using observations of the standard stars HR 6348 (KO III), HD 166780 (K4 III) 
and HD 173511 (K5 III). While all three standard stars were used for LL, only the first 
was used for SL. Where available, the mean flux from multiple data collection events 
(DCEs) for a single module/order /nod was computed, with the standard deviation 
from the mean providing a measure of the corresponding uncertainty. Spikes not elim- 
inated by replacement of bad pixel values were effectively removed by replacing any 
uncertainty larger than the average in the local wavelength range by the uncertainty 
from the nod position closer to the local average. Bonus-order spectra (corresponding 
to a small segment of the first-order spectrum obtained by the second-order slit) were 
averaged with spectra from the second and first orders over wavelengths they shared 
in common. 

The IRS spectrum for OGLE 051306.52-690946.4 (Fig. [32D shows typical features 
of carbonaceous dust - a featureless continuum at infrared wavelengths combined with 
a broad 11.3 /im SiC emission feature. There is also a 13.7 /im C2H2 feature. The 
strength of this feature ca nnot be reproduced b y a model photosphere alone, it is 



circumstellar in origin (see iMatsuura et al.ll2006l ) . We model both the circumstellar 



^For details, please refer to the IRS Data Handbook, available at 
http: / /ssc.spitzer.caltech.edu /IRS /dh/ dli32.pdf 

•^These modules are available in SPICE, the Spitzer IRS Custom Extraction package. 
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dust as well as the molecular gas in the extended atmosphere in §3.3[ 




Figure 3.2 The SED for OGLE 051306.52-690946.4 showing the mean OGLE I band 
flux (solid circle) and its range of variation (error bar), 2MASS (open triangles), 
IRSF (filled triangles) and SAGE (Epoch 1: open squares. Epoch 2: filled squares) 
photometry. The SAGE-Spec IRS spectrum (solid curve) is scaled down to fit the 5.8 
fim SAGE Epoch 1 photometry point. 

The optical, NIR and MIR photometry as well as MIR spectroscopy described 
above are incorporated into Fig. 13. 2[ The two-epoch photometry shows that there is 
no significant change in the shape of the MIR SED. We have taken advantage of this 
fact and scaled the SAGE-Spec spectrum down to the SAGE Epoch 1 5.8 /im flux to 
enable easy comparison to photometry. This scaling changes the flux predicted for 
the spectrum in the 5.8 fim band by ~5%. 
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3.3 Analysis 



3.3.1 2Dust Radiative Transfer Models 



The observed spectrum of an AGB star contains contributions from the central star 
as well a s the circumstellar dust shell. We model the AGB star+envelope system using 



2Dust flUeta k Meixner 



20031 ). a radiative transfer code for axisymmetric systems, 



with the simplifying assumption of a spherically symmetric dust shell. A comparison 
of the 2Dust output SED with the observed spectrum enables us to determine the 
AGB star mass-loss rate. 

As input, 2Dust needs information about the central star (effective temperature, 
size, distance and SED), the shell geometry (inner and outer radius, density varia- 
tion), and the dust grains (species, optical depth at a referei ice wavelength, grain 



size d istribution). We use the C-rich model photospheres of iGautschy-Loidl et al. 
( I2OO4J ) to represent the central star. To simplify the process, we choose a solar-mass, 
solar-metallicity model with an effective temperature of 3000 K and surface gravity 
\ogg = —0.4 and C/0 ra tio of 1.3 and place it at 50 kpc, the approximate distance 
to the LMC flFeastlligggf ). 

We assume a constant mass-loss rate and a constant outflow velocity Vcxp, leading 
to an inverse- square density distribution in the shell. Whi l e AGB outflow veloc ities 



are well-studied in the Galaxy (see, e.g., iLoup et al 



1993 



Olofsson et al 



19931 ). no 



measurements of outflow velocities exist for LMC carbon stars. While the lower 
metallicity of the LMC suggests a lower expansion velocity for 0-rich AGB outflows 
than for their Galactic counterparts, the outflow velocities of carbon stars are less 



sensitive to metallicity (jHabingj 119961 ). For simplicity, we assume a value of fexp=10 



km s~^ and ignore these complications. The ratio of the outer radius of the dust 
shell to the inner radius is kept fixed at 1000. The outer radius determines the 
total amount of mass in the shell, which in turn is related to the timescale of the 
mass loss. While this timescale is an important quantity, we are more interested in 
obtaining the AGB mass-loss rate, which is only weakly sensitive to changes in the 
outer radius as long as it is large enough to include contributions to the flux at the 
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longest wavelengths of interest (in our case, the 24 /im band) from the coldest dust 
in the shell. The mass-loss rate is, however, very sensitive to the value chosen for the 
inner radiu s. Observations of Galactic carbon stars {e.g., IRC -l-10216, 



1990 



19951) and results of radiative transfer modeling {e.g. 



Danchi et al. 



van Loon et al. 



1999 



Groenewegen et al.l Il998l ) suggest that the inner radius varies in the range ~2-20 



stellar radii, while simple estimates based on energ y balance sug gest that amorphous 
carbon dust should form within a few stellar radii (iHofnerl 120071 ). We vary the inner 
radius within these bounds until the desired SED shape is obtained. 

2Dust also requires us to specify the optical depth of the dust at a reference 
wavelength, which we choose to lie near the center of the SiC emission (11.3 yum). 
This optical depth is adjusted until a good fit to the SiC feature is produced. The 
spectrum of OGLE 051306.52—690946.4 has a strong 11.3 /im SiC feature, and a 
long-wavelength continuum. We model the dust around this star with amorphous 
carbon (optical constants from 



Zubko et al 



constants from 



19961) and silic on carbide (a-SiC, optical 
(120081) discuss the limitations of 



Pitman et al. 



using the a-SiC optical constants and provide newer, more reliable values, but they 
also require non-spherical grains. This is a complication we would like to avoid for 
the time being. 

The code requires a distribution of grain sizes from which it calculates cross- 
sections for an "average" grain. Our grain sizes are distributed according to the 
(IKim et al.lll994j ) (KMH) prescription, with sizes between 0.01 /im and ~1 /im (There 
is no strict upper bound on the KMH grain sizes; the latter number is an exponential 
scale height in the distribution.) Given this distribution of sizes as input, 2Dust then 
internally calculates absorption and scattering cross- se ct ions and asymmetry factors 
for an "average" spherical dust grain from Mie theory (jMielll908[) assuming isotropic 
scattering. We run the code in the Harrington averaging (IHarrington et al.l Il988l ) 



size chosen for dus t grains in radiative transfer models {e.g., 



is typically the 


single 


van Loon et al. 


1999; 



GroenewegenI |2006| ) . 



2Dust is now run for different values of the inner radius and 11 fim optical depth 
until the observed shape is reproduced. The fiux expected in the 5.8 /im band from 
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this model SED is then scaled to the 5.8 fim flux of the SAGE Epoch 1 data for OGLE 
051306.52—690946.4. The scaling changes the luminosity of the model photosphere 
as well. At constant effective temperature, this implies a change log (7. Such a change 
would in reality also change the shape of the model photosphere slightly, but this 
scaling is unavoidable given the small range of available models. The 2Dust input 
and output parameters are listed in Table 13.11 We will summarize the results of the 
RT modeling in §3.41 



3.3.2 Modeling The Circumstellar Gas 



Our model is also unable to reproduce some absorption features in the spectrum, 
despite the fact that our photospheric model takes into account the absorption due 
to molecules forming near the photosphere, leading us to believe that the strong 
absorption features seen in the data are a result of gas in the extended atmosphere 
and the circumstellar shell. We present a treatment of the circumstellar gas in this 
section. OGLE 051306.52—690946.4 has a few molecular signatures in its spectrum, 
such as C2H2 absorption at ~7 and 13.7 fim. The most prominent among these is 
the 13.7 fim absorption. We simulate the molecular band s of CO and C2H2 in OGLE 



051306.52-6 90946.4 using th e mode ls of 



Matsuura et al. 



published by 



Rothman et al. 



(120021 ) with recent line lists 



(120091). Assuming that the majority of the molecular 
absorption comes from the circumstellar shell, the excitation temperatures are chosen 
to be cooler than the temperature at the inner shell radius, i.e., Tex < ^in- Based on 
these considerations, we derive excitation temperatures of 1000 K for both species, 
with column densities of 10^^ (C2H2) and 3 x 10^^ (CO) cm~^ respectively. Fig. 13.31 
shows the locations of the spectral features for the parameters above. The resulting 
(gas+dust) fit to the spectrum is shown in Fig. 13.41 



54 



Table 3.1. 2Dust parameters for OGLE 051306.52-690946.4 



Central star'* 



Dust shell 



Dust grains 



Fit results 



Teft 3000 K 
logg -0.4 
C/0 1.3 



-Rin 
Rout 



6 i?* 

1000 Ru 



density profile p{r) 



10 km s 



Species 
Abundance 
r(11.3 /xm) 



AmCV SiC 
15% SiC 
0.25 

Size distribution KMH'^ 

(a,nin=0.01 ^m, 
ao=l fim, 7=3.5) 



L 5670 Lq 

Md 2.6x10-^ M0yr-i 

Tin 1150 K 



'^Photosphere model from iGautschv-Loidl et al.l (|2004f ). 
Amorphous carbon grains, p — 1.8 g cm~"^, optical constants from Zubko et al. (1996 



■^a-SiC grains, p = 3.22 g cm optical constants from iPegourie ( 1988f ). 
'^Size distribution from lKim et al. (1994), n(a) ~ a"''' exp (— a/ao). 




Figure 3.3 The molecular features used to model the gas in OGLE 
051306.52—690946.4. The C2H2 and CO traces are shifted down for clarity. 



3.4 Results 

We find that the best-fit model SED is produced for an inner radius six times 
the stellar radius and an optical depth of 0.25 at 11.3 fim. The luminosity of the 
model photosphere is scaled up by a factor of 1.29 to 5670 Lq to fit the observations. 
As mentioned in §3.3.11 a higher luminosity at constant effective temperature can 
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Figure 3.4 The results of modeling the gas and dust around OGLE 
051306.52-690946.4. The SAGE Epoch 1 (red dots) and 2 (black dots) photome- 
try as well as SAGE-Spec spectrum (black curve, scaled to Epoch 1 photometry) 
shown with the best-fit dust RT model, scaled to fit Epoch 1 and Epoch 2 photom- 
etry (cyan curves). The model SED is convolved with the molecular features shown 
in Fig. 13.31 to show the effect of adding 02^.2 only (pink), CO only (green) and the 
combined effect (blue). 

only be obtained by changing log (7, which will affect the photospheric spectrum. We 



disregard this complication 



L (see, e.g.. iGroenewegen 



or thi s study. Assuming that the mass-loss rate scales as 
20061 ). we find the dust mass-loss rate to be 2.6 x 10~^ 



M0yr~^. Since the shell geometry is defined in terms of the stellar radius, we do not 
expect the dust temperature to change as a result of scaling. The temperature Tin of 
the dust at the distance corresponding to the inner radius -Rin of the dust shell is 1160 
K, which is close to the condensation temperature of a-SiC. We also fit the 2Dust 
output to the second epoch of SAGE photometry with good agreement. In this case, 
we obtain a luminosity of 5330 LQand a corresponding mass-loss rate of 2.2 x 10^^ 
Moyr-i. 



57 



Fig. 13.41 shows the 2Dust output SED (cyan curve) superimposed on the SAGE- 
Spec spectrum and Epoch 1 photometry (red dots). While the fit is acceptable overall, 
we need to incorporate the effects of circumstellar gas to fit the 13.7 /im feature. The 
figure also shows a convolution of the 2Dust output SED with the molecular features 
(C2H2 only: pink; CO only: green; C2H2 + CO: blue), which is in excellent agreement 
with the SAGE-Spec spectrum. 

We obtain an estimate for the gas mass-loss rate from our results for C2H2 (the 
CO band is at the edge of our spectrum). Depending on the location of the C2H2 line 
forming region, the mass-loss rate ranges between 1 x 10""^ - 2 x 10~ ^ M^vr"^ Here, 
we h ave adopted a C2H2 abundance of 10~^ compared to H2 (see 



Matsuura et al. 



20061 ). but some of the C2H2 will have condensed out as dust. 



3.5 Discussion 



OGLE 051306.52-690946.4 is a 



1 yr (llta et al. 



2004 



Groenewegen 



ong-p eriod variable (LPV) with a period of about 



20041 ). The two epochs of SAGE observations 



have thus captured the star at different stages in its pulsation. While 2MASS and 
IRSF observations in the NIR also offer multiple-epoch data, these were not taken in 
conjunction with SAGE data. The amplit ude of variation is expected to be higher 



in the NIR magnitudes than in the optical (IGroenewegen 



2004. estimated the I-band 



amplitude to be about 0.8 mag). The dust mass-loss rate predicte d by our 2Dust 



mode 



van Loon et al. 



2.6 X 10~^ MQyr^^, is on the low end of the rates derived by 
( ll999l ) for the brightest, most obscured carbon stars in the LMC. The derived dust 
MLR can be converted to gas MLR of 5.2 x 10"'' M0yr~^, assuming a gas:dust ratio 
for C-rich AGB stars of 20 0. This is comparab l e to t he ga s MLR estimated from the 
period vs MLR relations of IVassiliadis fc WoodI (1l993l ) and IGroenewegen et al.l (1l998l ) 
( 10""^ and 7.6 x 10~^ Mf^yr~^respectivel y, assuming a pe r iod of 358 d) and the MLR 
predicted by the color vs MLR relation of 



Matsuura et al. 



mm (4-6x10-^ Mqjt 



-1 



using Kg- [8.0 
relation from 



-0.395 mag for OGL 



Srinivasan et al. 



051306.52-690946.4). The 8 fim excess - MLR 
(j2009|) gives 2 x 10"^ MqJt-^ and 1.26 x 10"^ MqJt-^ 
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for the Epoch 1 and 2 observations respectively. Despite the fact that the Kg and 8 
fj.m data were obtained for different phases, the MLR predicted by the Ks-[8.0] color 
agrees well with the MLR determined from other methods. 

The long-term goal of this study is to use 2Dust modeling of carbon star circum- 
stellar shells to fit the photometry of the ~7000 C-rich AGB candidates in the SAGE 
survey. In this section, we weight the simplifying assumptions used in our modeling of 
OGLE 051306.52—690946.4 against this future goal. Our combined gas+dust model 
improves on the spectrum predicted by 2Dust modeling of the dust shell alone. There 
are, however, many issues associated with our simple assumptions in this work. For 
instance, our model requires a significant amount of SiC (15%). Furt hermore, the 



value of Tin is very close to the condensation temperature of a-SiC. ( iSpeck et al. 



20091 ) have shown that the use of graphite over amorphous carbon can result in lower 



amounts of SiC needed to fit the observed 11 /im feature, as well as lower tempera- 
tures for the dust. 



From Fig. 13. 4[ it is clear that our AmC+SiC dust composition reproduces the 11.3 
/im feature and the continuum in overall observed spectrum well. However, the SiC 
feature is observed to be broader than the model prediction; this is in p art a problem 



with the 



Pegourid (119881 ) optical constants for a-SiC, as discussed by 



Pitman et al. 



(120081), who suggest the usage of /3-SiC optical constants for non-spherical grains. 
The observed spectrum also looks broader than those of known Galactic and LMC 
carbon-rich spectra. A spike noise may be increasing the shape artificially. For the 
purpose of modeling the SAGE photometry, however, the details of the SiC feature 
are irrelevant as the SiC feature does not fall in any of the Spitzer bands. 

The model also systematically overestimates the flux beyond ~ 20/im. Although 
the quality of data in the long-wavelength range is poor, this disagreement may also 
be due to the simplifying assumption s in our model. Our assumption of spherical 



grains is also an idealization, and the 



Kim et al 



(1l994l ) prescription may not repre- 



sent the true distribution of dust grain sizes in the cir cumstellar shells o f AGB stars, 
but these details do not affect the overall properties (jSpeck et al.l 120091 ). Our main 
interest is to derive dust composition and mass-loss properties for the entire sample 
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of SAGE C-rich AGB star candidates using a baseline dust description, and for these 
purposes the assumption of spherical grains is sufficient. 



3.6 Summary and conclusions 

We have modelled the circumstellar shell a carbon-rich asymptotic giant branch 
star using the radiative transfer code 2Dust to fit SAGE photometry and SAGE-Spec 
spectroscopy for the source OGLE 051306.52—690946.4. We find a dust mass-loss 
rate of 2.6x10"^ Mgyr"^ for the source, corresponding to an inner radius for the dust 
shell of six times the stellar radius. The optical depth at 11.3 /xm is found to be 
0.25 and the dust temperature at the inner radius of the shell is 1160 K. We model 
the molecular features in the circumstellar shell with slab models of C2H2 and CO. 
The gas mass-loss rate estimated from the column density of C2H2 is in the range 
1 X 10~^ — 2 X 10~^ M0yr~^and it depends heavily on the location of the line-forming 
region. The model SED when combined with C2H2 and CO models shows excellent 
agreement with the observed spectrum and photometry. We are successful in our aim 
of defining a baseline set of dust properties for use in our model grids, which will be 
the focus of the following chapter. 
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Chapter 4 

A Grid of Carbon Star Models 



4.1 Introduction 



Asymptotic giant branch (AGB) stars are among the bright est IR objects in 



star-forming ga" 



Srinivasan et al. 



axies. The ~ 45 000 AGB candidates (see, e.g., 



Blum et al. 



2006 



20091 ) inject a significant fraction of their mass (30- 80%) in the form 



of ga s and dust into the LMC ISM at a total rate of ~ 10 ^ M0yr ^ ( ISrinivasan et al. 



20091 ). thereby increasing its metallicity and driving the chemical evolution of the 
galaxy. AGB stars are thus essential ingredients in modeling galaxy evolution. 

Despite being outnumbered by the 0-rich AGB stars, carbon stars are probably 
responsible for the majority o f the AGB mass ejectio n in the LMC because most 



extreme AGB stars are C-rich. 



Srinivasan et al 



(120091 ) estimated the total mass loss 



return from the near- IR identified C-rich AGB stars to be 4.8 x 10"^ Meyr"^ and 
the extreme AGB stars to be (4.7 — 11.8) x 10^^ MQyr^-*^. However, their estimate 
was an empirical one, based on excess mid-IR fiux over the photosphere and using an 
excess versus mass-loss rate relation derived from data for a handful of C-rich and 
extreme AGB stars. Fits to spectral energy distributions for the AGB stars using 
radiative transfer model grids provides a more direct and quantitative measure of the 
mass-loss rate return by these stars. 



Recent work has resulted in many model grids for carbon stars: iMarigo fc Girardi 



( 120071 ) presented models of synthetic thermally pulsing (TP) AGB evolution and de- 
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rived theoretical isochrones from these models in 



Marigo et al. 



( 20091) improved u pon the dynamical wind models of 



(2008). 



Mattsson et al. 



Mattsson et al. 



Hofner et al 



20071 ) and 



(120031 ) a nd computed 900 carbon star models to study the ir dependence 



on ste llar parameters. iGautschy-Loidl et al.l (120041 ) and, more recently, lAringer et al. 



(120091 ) constructed a grid of hydrostatic photosphere models for C-rich AGB stars. 
In this chapter, I introd uce a grid of carbon-s tar models that calculates the effect 



of circumstellar dust on the 



Aringer et al.l (20091) mo del photospheres using the dust 



radiative transfer code 2Dust ( lUeta fc Meixnerll2003l ). The models span a wide range 
of stellar and dust-related parameters. I also present synthetic photometry for these 
models in the 2MASS, Spitzer, AKARI and WISE bandpasses. In §4.21 1 discuss the 
various input parameters that went into the construction of the grid. I briefly describe 
the details of the method for computing the models in §4.3[ and discuss the resulting 
model grid and synthetic photometry in §4.41 I present my conclusions in §4.51 



4.2 Input to the radiative transfer code 



4.2.1 Model photospheres 



As input to 2Dust, I utilize the lAringer et al.l (120091 ) photosphere models, which 
contain 132 models with meta llicity oi Z = 0.33, close to that of the LMC {Z ~ 
0.3 — 0.5 Zq, IWesterlundl 119971 . p. 234). Their LMC grid was calculated completely 



for 2 Mq models, with a subset of 1 Mq models for Tefj=2600 K. The choice of C/0 
ratios (1.4, 2.0 and 5.0) reflected the fact that this ratio can increase above unity in 
low-metallicity environments such as the LMC after a single thir d dredge-up event. 
C/0 ratio for more evolved s tars c an rise significantly above unity (iMarigo fc Girardi 
20071 ). While lAringer et al.l (120091 ) considered a temperature grid of 2600 K to 4000 
K in increments of 100 K, we ignored the models with T^s values higher than 3500 
K as these temperatures correspond to the post-AGB phase. This restriction on Tgfj 
reduces the number of available photosp here models to 120. C arbon stars may have 
effective temperatures as low as 2400 K (IMattsson et al.l 120091 ). However, the lower 
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temperature limit of 2600 K for the lAringer et al. 



2009) 



coverage for the LMC C-rich sample (jGroenewegen et al. 



mode ls provides adequate 



200% . 



The surface gravity ranged from log (7=0 down to log(7=-l in steps of 0.1, but 
not all values were available for each T^s value because of convergence issues. For 
Teff=2600 K, the only \ogg=0 was considered. The corresponding luminosity is about 
1000 L(7), which is significa ntly less than that for typical carbon stars (see, e.g., Fig. 



4 in 



Srinivasan et al. 



20091 ) . The lower limit of log g d epends on Tpff. It was ch osen 



to cover the region of AGB carbon stars predicted by iMarigo fc Girardil (120071 ) and 
Marigo et al.l (120081 ). For the range of model param eters corresponding to each tem- 
perature range, see Table 1 in lAringer et al.l (120091 ). The luminosities of the selected 
models range from 10^ LQto 3 x 10^ L©. The most luminous optically visible carbon 
stars are at luminosities slightly higher than this upper limit, and the extreme carbon 
stars are perhaps as luminous as ~ 10^ L©. The lack of more luminous model photo- 
spheres means that in order to model the most evolved, highly enshrouded stars, we 
will have to scale the luminosities of the available models. For the optically thick dust 
shells, the scaling of luminosities should provide adequate matches to the available 
data which measures the total luminosity of the star+envelope system. 



4.2.2 Circumstellar envelope 

In chapters [2] and [3], I probed the parameter space for modeling by first estimating 
color temperatures and optical depths of the dust, followed by the radiative transfer 
modeling of a star for which photometry and spectroscopic data were available. In this 
section, I will use these results to constrain the parameter space for the circumstellar 
envelope. 



4.2.2.1 Shell geometry 

With the simplification of spherical symmetry, we are left to specify the inner and 
outer radii of the shell, as well as the density structure in the shell and the expansion 
velocity of the dust. The outer radius is kept fixed at a thousand times the inner 
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radius. As mentioned in chapter [3l this is acceptable as long as we do not use it 
to estimate the timescale of mass loss. Four values are chosen for the inner radius, 
at 1.5, 4.5, 7 and 12 times the radius of the star. This range of ra dii samples th e 
theoretically expected range of dust production in carbon star winds (jnofner boOT ). 



The expansion velocity is fixed at 10 km s ^. 



wind should in general be luminosity-dependent (IGail fc Sedlmayr 



he velocity of a radiat ion-driven 



198611 and this 



dependen ce can be accou nted for by an appropriate scaling relation with luminosity 



[see, e.g., 



van Loon 



20071 ). However, in this first attempt at generating a model grid. 



we will ignore this complication. 

We also choose an inverse- square density drop-off for the shell. This implies a 



constant mass-loss rate, which may be unrealistic over the cour s e of AGB evo 



1987 



Zijlstra et al 



ution 



1992 



or during rapid th ermal pulsing. Some authors (iBedijn 
Groenewegenlll993l ) have considered time- varying mass- loss rates to account for such 
situations. Since our main interest is fitting a large set of photometric data in or- 
der to derive global mass injection rates, this is beyond the scope of our present 
study. 2Dust is capable of accepting user-defined density functions, however, and an 
investigation of density profiles different from can be undertaken if desired. 



4.2.2.2 Dust properties 

As discusse d in chapter [3l w e choose a mixture of amorphous carbon (optical 



constants from 



Zubko et al 



and silicon carbide (a-SiC, optical cons tants from 



Pegourielll988l) . with 15% SiC. The grain sizes are distributed according to 



Kim et al. 



( 1l994j . KMH) with a minimum grain size of 0.01 /im and a "maximum" size of 1 
/im. We vary the optical depth at the center of the SiC feature at 11.3 /im from 
Til. 3 ^lm = 10^^ to 1 in equally spaced logarithmic steps with five values per decade. 
Furthermore, we also consider four th 3 ^m=2, Tefi=2600 K models. 

This range covers most of the optically visible carbon-rich AGB stars, as well as 
the lower luminosity extreme AGB stars. We will probe higher optical depths in future 
studies in order to model the most extreme carbon stars. The dust temperature at the 
inner radius is another parameter that we constrain in our modeling, we only accepted 
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Table 4.1. Input parameters for model grid 



Parameter 



Value or Range of values 



Luminosity 
Inner shell radius 
Outer shell radius 
Density profile 
Expansion velocity 
Dust species 



- fO^ to 3 X lO'^ Lq 
1.5, 4.5, 7, 12 
1000 



10 km s ^ 
85% AmCi- 



15% SiC2 



Optical depth at 11.3 /im 10 to 1 (five values per decade)"^ 



^Optical constants from lZubko et al.l ( 19961 ) 
Optical constants from iPegourid (|l988l ) 
^Additionally, four rcff=2600 K, t==2 models. 



models that resulted in a value of T(i?in) lesser than or equal to the condensation 
temperature of a-SiC is ~1800-2000 



With the above constraints on the 



Aringer et al.l (120091 ) photosphere parameters 



and the various dust parameters, about 7750 models are available to run through 
2Dust. The parameter ranges for these models are summarized in Table I^TTl Placing 
an added constraint on the dust temperature brings the number of models down to 
about 7100, as discussed in the following section. 



4.3 Computational details 

The models were run on the Royal Linux computing cluster at the Space Telescope 
Science Institute. The cluster consists of 23 computing nodes, each with 2 x 2 x 2.4 
GHz processors, 8 GB RAM, 160 GB local disk space and 17 TB of total disk space. 
Royal enables parallel processing at high computing speeds, and also automates the 
job submission and scheduling. 2Dust was therefore run in the non-interactive mode. 
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with an initial run to ensure that the number of radial grid points was sufficient. 
The 120 photospheres were split into at most 19 photospheres per job submission, 
totahng up to ~ 1200 individual models submitted simultaneously. The time taken 
for successful convergence ranged from ~ 10 min for the optically thin models to a 
few hours for the models with moderate optical depth. The Royal cluster had an 8 
hour time limit on individual submissions, which was insufficient for optical depths 
higher than unity. A smaller subset of r = 2, Tefj=2600 K models with log(7=-0.9, 
-0.8, -0.6 and -0.5 were executed successfully for two values of inner radius (-Rin/-R* 
= 7 and 12) on a Mac Pro Desktop with a 2 x 2.8 GHz Quad-Core Intel Xeon CPU 
and 10 GB RAM. In general, the code required an increase in the number of radial 
grid points to deal with higher optical depths, which eventually resulted in memory 
allocation errors. At present, therefore, our models are reliable up to r = 2. Once the 
models converged, only those runs that produced dust temperatures lesser than 2000 
K were considered as successful. This resulted in the rejection of some i?in = 1.5-R,, 
models, which for hotter stars naturally produced warmer dust temperatures. 

I then produced synthetic photometry for the remaining models in the optical 
UBV10, 2MAS^ and Sj>zfee?§ bands in order to compare t he results with SAGE 
data. Synt hetic photometry wa s also derived for the AKARI (IMurakami et al.l 120071 ) 



and WISE (IWright et al. 



200J) passbandf^. 



^The MCPS magnitudes were placed on the Johnson-Kron-Cousins UBVI system. The de- 
tector quantum efficiency curve was obtained from the Las Campanas Observatory website, 
http://www.lco.cl/lco/index.html. Filter profiles for the Johnson U , Harr is B, V, and Cousins I 
filters were obtained from the references in Table 9 of iFukugita et al" I (|l995l) . 

^The 2MASS filter relative spectral responses (RSRs) derived by ICohen et all (|2003( ) 
were obtained from the 2MASS All-Sky Data Release Explanatory Supplement, available at 
http://www.ipac.caltech.edu/2mass/re leases/allsky / doc/s ec64a.html. 

■^The IRAC RSRs are plotted in iFazio et al.l ( 2004 ). and were obtained from the Spitzer 
Science Center IRAC p ages at http://ssc.spitzer.caltech.edu/irac/spectralresponse.html; The 
MIPS ( Rieke et al.l 120041 ) RSRs were obtained from the Spitzer Science Center MIPS pages at 
http://ssc.spitzer.caltech.edu/mips/spectralresponse.html 

''Filter response curves available at 
http://www.ir.isas.jaxa.jp/ASTRO-F/Observation/RSRF/IRCFAD/index.html 
and http:/ /www. astro. ucla.edu/~wright/WISE/passbands. html respectively 
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4.4 Results and discussion 



In this section, I demonstrate the success of our model grid in covering the range of 
luminosities and mass-loss rates expected for carbon stars. I use synthetic photometry 
to demonstrate the coverage of our model grid in color-color and color-magnitude 
diagrams (CMDs) and compare it to SAGE data for AGB star candidates. I also 
provide color-color diagrams for the AKARI and WISE bands, and I revisit the model 
for OGLE 051306.52—690946.4 and compare my results in Chapter [3] to the best-fit 
model produced by using the model grid generated in this chapter. Finally, I show 
preliminary results of fitting the model grid to SAGE C-rich AGB candidates. 

4.4.1 Dust mass- loss rates 

Fig. 14.11 shows the range of luminosities and mass-loss rates explored by the 
current model grid. The available luminosities span the range of C-rich AGB star 
luminosities observed in the SAGE survey (see Fig. 12.41 in Chapter [2]). The highest 
mass-loss rates are ~ 10"'' Moyr^^, corresponding to a total rate of 2 x 10~^ MqYi~^ 
(assuming a gas:dust ratio for carbon dust of 200). This covers the range of mass-loss 
rates expected for the low- and moderately-obscured carbon stars. More models for 
higher optical depths will be required to reproduce the mass-loss rates of the extreme 
AGB stars. 

4.4.2 SAGE synthetic photometry 

The models are superimposed on SAGE C-rich and extreme AGB candidates in 
Fig. 14. 2[ The pink line shows how the Tcfj=2600 K, log (7 = —0.9, Rni = 7-Rstar model 
would evolve with increasing optical depth, assuming that the rest of the parameters 
stay fixed. Fig. 14.31 shows a similar variation. The IRAC-IRAC color-color diagram 
(Fig. 14. 4p shows that we underestimate the [5.8]-[8.0] colors of the optically thin 
C-rich AGB stars or lack coverage of bluer C-rich AGB stars. This could be due to 
the presence of small grains in our models - the minimum grain size was set to 0.01 
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Figure 4.1 The range of values covered in luminosity- mass-loss rate space by the model 
grid. The four r=2 models and the r=6, 8 and 10 models can be clearly seen. The 
dashed line represents the classical AGB luminosity limit. 



/im, whereas meteoritic evidence points to the fac t that ~70% of SiC grains have 



sizes in the range 0.3-0.7 fim (see ISpeck et al. 



20091 . and references therein). We will 



investigate the effect of grain sizes on such deviations in detail in follow-up studies. 
The figures show that the models reproduce the colors of the low- and moderately- 
obscured carbon stars to a great extent, despite the simplifying assumptions that 
went into creating the grid. The grid allows us to readily model the mass-loss rates 
of carbon stars of low and moderate optical depths. 

Knowledge of the distance to the LMC allows us to also construct color-magnitude 
diagrams (CMDs) for our models. Fig. K5[ WB and lO show the 2MASS, IRAC and 
IRAC-MIPS CMDs comparing the models to SAGE C-rich and extreme AGB can- 
didates. Once again, we see that the models reproduce the data well. The [8.0] vs 
[3.6]-[8.0] CMD shows that beyond a certain optical depth, the 8 /xm magnitude is 
constant for a large range of colors. This corresponds to the spectral energy distribu- 
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Figure 4.2 The Ks-[8.0] vs Ks-[3.6] color-color diagram showing the model grid (black 
dots) superimposed on the C-rich (red dots) and extreme (green dots) AGB candi- 
dates. The pink line shows the track for a single model (Teff=2600 K, logg = —0.9, 
i?in = 7i?star) with increasing optical depth. 

tion (SED) peaking in the 8 fim band. The 24 fim flux for the single track considered 
in Fig. 14.71 undergoes a sudden drastic increase at almost constant [8.0]-[24] color. 
This probably corresponds to the onset of appreciable mass-loss. 



4.4.3 Synthetic photometry for other surveys 

In this section, I provide color-color diagrams for the A KARI and WISE passbands 
generated from synthetic photometry of our model grid. Ilta et al.l (120081 ) presented 
the AKARI survey of the LMC and discussed color-color and color-magnitude dia- 
grams. The AKARI survey is important for AGB stars because the Sll band directly 
samples the 11.3 fim SiC feature in C-rich AGB stars as well as the ~ 10 /im sil- 
icate emission in 0-rich AGB s tars. The increa se in feature strength, therefore, is 



reflected in the [S11]-[L15] color. 



ItaetaL 



(120081 ) found by generating synthetic pho- 
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K-[3.6] 



Figure 4.3 The Ks-[24] vs Ks-[3.6] color-color diagram showing the model grid (black 
dots) superimposed on the C-rich (red dots) and extreme (green dots) AGB candi- 
dates. The pink line shows the track for a single model (Teff=2600 K, logg = —0.9, 
i?in = 7i?star) with increasing optical depth. 



Kraemer et al. 



(12002) that 



tometry in the AKARI bands for the ISO-SWS spectra of 
the change in the SiC feature strength manifests itself as first a blueward shift of the 
[S11]-[L15] color followed by a reddening. Fig. 14.81 shows the iKraemer et al.l (120021 ) 



sources classified as stars with carbonaceous molecular features (NC), stars showing 
carbon dust (CE/CR), and carbon-rich post- AGB objects (CN/CT) plotted on the 
grid of models. The model grid do es not show the same blueward-redward shift as 

(I2OO2I ) spectra. We note that in 



seen in the colors derived from the 



Kraemer et al. 



future studies, LMC data from AKARI will perhaps help constrain the role of SiC 
in our models better than SAGE data because of the direct detection of the 11.3 fim 
feature. The [S11]-[L15] color is also a tracer of mass-loss. 

Fig. mi shows the ratio of the model fluxes in the WISE W3 (12 /xm) and W4 (22 
fj,m) bands as a function of the corresponding ratios in the Wl (3.4 /xm) and W2 (4.6 
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[3.6]-[4.5] 



Figure 4.4 The [5.8]-[8.0] vs [3.6]-[4.5] color-color diagram showing the model grid 
(black dots) superimposed on the C-rich (red dots) and extreme (green dots) AGB 
candidates. The pink line shows the track for a single model (Teff=2600 K, logg = 
—0.9, Rin = 7i?star) with increasing optical depth. 

yum) bands. The WISE W3 band can directly detect the 11.3 /xm SiC feature and 
is therefore an important band for AGB studies. Similar to the AKARI Sll band, 
when WISE data is available, information from the W3 filter will help improve the 
treatment of SiC dust in our model grid. 

From these diagrams, it is clear that overall, our models provide good coverage 
over a wide range of parameters. The phenomenon of variability can also be studied 
using our model grids. The availability of information about the SiC feature in the 
AKARI and WISE surveys will help constrain the behavior of the dust shell over the 
course of the pulsation, which will assist us improving our model prescriptions. 
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Figure 4.5 The Kg vs J-Kg CMD showing the model grid (black dots) superimposed 
on the C-rich (red dots) and extreme (green dots) AGB candidates. The pink line 
shows the track for a single model (Teff=2600 K, \ogg = —0.9, -Rin = 7-Rstar) with 
increasing optical depth. 



In order to calculate mass-loss rates for all the SAGE AGB star candidates, we 
need to fit their photometry to that synthesized from the model grid. In order to 
verify this, I fit all four IRAC magnitudes and the MIPS 24 fim magnitude of the 
SAGE data to the models as a first attempt at obtaining best-fits. I minimized 
for each source where 



here, and Srrii are the magnitude and error of the AGB candidate in the ith 
wavelength band, mmod,i are the corresponding magnitudes for the models, and the 
sum is over all four IRAC bands and the MIPS 24 band. In this section, I show some 
results from this simple fitting technique in order to demonstrate the utility of the 



4.4.4 SED fits 




(4.1) 
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[3.6]-[8.0] 



Figure 4.6 The [8.0] vs [3.6]-[8.0] CMD showing the model grid (black dots) superim- 
posed on the C-rich (red dots) and extreme (green dots) AGB candidates. The pink 
line shows the track for a single model (Teff=2600 K, \ogg = —0.9, Rin = 7i?star) with 
increasing optical depth. 

model grid. 



4.4.4.1 Revisiting the model for OGLE 051306.52-690946.4 



In Chapter [31 1 used the 



Gautschy-Loidl et al. 



(120041 ) model photospheres to calcu- 



late radiative transfer models for the variable carbon star OGLE 051306.52—690946.4. 
In this chapter, I cor npare the pred i ctions of this model with those obtained from our 
model grid using the lAringer et al.l (l2009l ) photospheres. The results are summarized 
in Table 14.21 Th e model grid is able to predi ct a mass-loss rate very similar to that 
predicted by the iGautschy-Loidl et al.l (120041 ) models, and it also predicts a similar 
temperature for the dust at the inner radius. Both best-fit models correspond to the 
same Tegbut different surface gravities. This latter effect is pe rhaps due to the fact 



that there are a limited number of log g values available for the 



Gautschy-Loidl et al. 
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12 



[8.0]-[24] 



Figure 4.7 The [24] vs [8.0]-[24] CMD showing the model grid (black dots) superim- 
posed on the C-rich (red dots) and extreme (green dots) AGB candidates. The pink 
line shows the track for a single model (Teff=2600 K, \ogg = —0.9, Ri^ = 7i?star) with 
increasing optical depth. 



(I2OO4J ) models. The model in Chapter [3] was constructed with due attention paid to 
the spectrum of OGLE 051306.52—690946.4. Despite simply fitting to the photome- 
try of the source in this chapter, we have obtained a similar result for the mass-loss 
rate of the source. This demonstrates the effectiveness of the model grid. 



4.5 Conclusions 

I constructed a grid of carbon star models using the radiative transfer code 2Dust 
and for a range of various stellar and dust shell parameters. At present, the model 
grid only covers optical depths at 11.3 /im up to 2, so it can not yet be used to 
model the most obscured stars. In the future the coverage will be extended to higher 
optical depths. We will also address the effects of grain chemistry, sizes and shapes 
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Figure 4.8 The AKARI [N3]-[S11] color vs [S11]-[L15] color, comparing the cov- 
erage of the rnqdel g rid to the locations of spectroscopically classified stars from 



Kraemer et al.l (120021 ) showing carbon-rich features (red crosses: NC, green dots: 
CE/CR, pink dots: CN/CT). The pink line shows the track for a single model 
(Tefj=2600 K, \ogg = —0.9, i?in = 7-Rstar) with increasing optical depth. 

and investigate the dependence of the models on the metallicity of the environment. 

I provide synthetic photometry in the optical, 2MASS and Spitzer bands as well 
as in the AKARI and WISE bandpasses. The latter are chosen for the purpose of 
directly studying the variation of model parameters with the observed strength 11.3 
fim SiC feature. When AKARI and WISE data becomes available, this model grid 
will prove to be very helpful in determining the dust properties of and the mass-loss 
return from the evolved stars in these data sets. One of the aims of this thesis is to 
devise a general-use fitter that can be used for large sets of photometric data. 

A concurrent effort to model 0-rich AGB and red supergiant stars is being un- 
dertaken by Sargent et al. (in preparation). When combined with the C-rich models 
presented in this thesis, these models will enable the assessment of mass loss return 
from galaxy-wide point source catalogs from projects such as SAGE. While most 
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Figure 4.9 The ratio of the fluxes in the W3 and W4 bands as a function of the flux 
ratio in the Wl and W2 bands. The pink line shows the track for a single model 
(Tefr=2600 K, \ogg = —0.9, Rin = 7i?star) with increasing optical depth. 

extreme AGB stars are probably carbon stars, some are very bright OH/IR stars. 
As of yet, it is not possible with photometry alone to distinguish between these two 
chemical types. The two grids will be used in conjunction to define the separation of 
C-rich and 0-rich sources in this extreme AGB star regime. 
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Table 4.2. Comparison of 2Dust fits to OGLE 051306.52-690946.4 



Quantity^ Chapter [3] This chapter 



Photosphere model Gautschv-Loidl et aL (2004) Aringer et aL (2009) 





3000 


3000 


L (Lo) 


5670 


6166 (C) 


logg 


-0.51 (Ltd.) 


-0.2 


M (Mq) 


1.0 (Ltd.) 


2.0 (Ltd.) 


-Rin (R*) 


6 


4.5 (C) 


Til. 3 /jm 


0.25 


0.2 (C) 


M (Moyr-i) 


2.57x10^^ 


2.16x10-^ 


r(i?in) (K) 


1160 


1140 



"'^We denote by (Ltd.) or (C) the model parameter values that were 
not constrained well owing to limited availability of models or coarseness 
of grid, respectively. 
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Chapter 5 



Concluding Remarks 



The main aim of this thesis was to be able to quantify the dust mass-loss return 
from asymptotic giant branch (AGB) stars to the Large Magellanic Cloud (LMC). 
This determination requires that we model the circumstellar dust shells around these 
stars individually, which is not easy considering the sheer number of AGB stars in 
the LMC. One way to achieve this is to generate a grid of models spanning a range 
of stellar and dust shell parameters typical of stars on the AGB which can then be 
fit to observed photometry in order to derive mass-loss rates. I have presented such 
a model grid in this thesis. These models are in progress of being extended as well as 
improved based on comparisons to data^ 

Using data from the SAGE survey (IMeixner et al.ll2006l ). I first identified 0-rich, 
C-rich and extreme AGB star candidates based on their locations in various near- and 
mid-infrared (IR) c olor- magnitu de diagrams (CMDs). Knowledge of the distance to 



Feast 



19991 ) enables us to also estimate the luminosities of these 



the LMC (-50 kpc, 

AGB candidates, allowing a study of the luminosity function of the sample. A future 
paper will address this topic in detail. SAGE data was also able, for the first time, 
to identify a faint, red population of 0-rich AGB stars. The theoretical isochrones 



of 



Marigo et al.l (120081 ) explain that the bright population corresponds to massive 



sources that undergo hot-bottom burning and stay 0-rich, while stars in the faint 
population (which consists of about 80% of the 0-rich stars) undergo third dredge-up 
and eventually turn C-rich. I compare the spectral energy distributions (SEDs) of the 
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SAGE AGB candidates with model photospheres in order to calculate their mid-IR 
excesses, and use these excesses in turn to estimate the color temperature and optical 
depth of their dust shells. I find that the IR excess increases with luminosity for all 
three types of sources in all four IRAC bands as well as the MIPS24 band. The color 
temperature shows a decrease with increasing excess, and the optical depth shows an 
increasing trend. I conclude that the IR excess is a good proxy for the mass-loss rates 
of AGB stars. I quantify this relation by deriving a fit to the excess-mass-loss rate 
relation for a handful of well-studied AGB stars and using this relation to calculate 
mass-loss rates for all the sources in order to estimate the total dust mass loss return. 
I find that the dust injection rate from AGB stars is about (5.9-13) x 10^^ MQyr^^. 

In order to constrain the dust properties that will be incorporated into the carbon 
star model grid, I perform 2Dust radiative transfer modeling of the dust around 
OGLE 051306.52—690946.4, a variable carbon star for which both photometric and 
spectroscopic data is available through the SAGE and SAGE-Spec (Kemper et al. 
2009, in preparation) studies, respectively. The source has a prominent 11.3 fim SiC 
emission feature, which I use to constrain the ratio of SiC to amorphous carbon used in 
the modeling. I find that with 15% of SiC, I find a good fit to the spectrum of OGLE 
051306.52—690946.4, resulting in a dust mass-loss rate of 2.6 x 10"^ Moyr^^. We 



also model the CO and t 



models (jMatsuura et al. 



le 13 .7 fim C2ll2 feature in the circumstellar shell with slab 



2OO2I ) . From the strength of the C2H2 feature, we estimate 



a gas mass-loss rate in the range 1 x 10^ ''-2 x 10~^ M0yr~^, which depends heavily 
on the location of the line-forming region. 

Finally, I present a grid of carbon star models constructed with the radiative 
transfer code 2Dust for a range of various stellar and dust shell parameters. I provide 
synthetic photometry in the 2MASS and Spitzer bands as well as the filters used in 
the AKARI and WISE surveys. The color-color and color- magnitude diagrams show 
that despite our simplifying assumptions we are able to reproduce the general features 
of the SAGE photometry well. A preliminary fitting scheme is then used to find the 
best-fit model to the photometry of OGLE 051306.52-690946.4, and I find that this 
best-fit model agrees well with that predicted in the Chapter [3l At present, the model 
grid reaches to an optical depth of 2 at 11.3 fim, so it can not be used to fit the most 
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obscured extreme AGB stars. In the future this coverage will be extended to higher 
optical depths. Along with a companion grid of 0-rich AGB models (Sargent et al., 
in preparation), these models should help distinguish 0-rich and C-rich features in 
extreme AGB stars. We will also consider the effect of changing the grain chemistry, 
grain size and non-spherical shape. The AKARl and WISE bandpasses contain filters 
that can directly detect the 11.3 /xm SiC emission, making them important for AGB 
star studies. When these data become available, we will be able to better constrain 
the dust properties used in the model grid. The model grid will also incorporate 
a general-use fitter that can be used for large sets of photometric data, enabhng 
accurate determinations of mass loss return from photometric data of entire AGB 
populations in galaxies. 
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